Recoil Polarization Measurements for Neutral Pion Electroproduction at 
Q 2 = 1 (GeV/c) 2 Near the Delta Resonance 
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We measured angular distributions of differential cross section, beam analyzing power, and recoil 
polarization for neutral pion electroproduction at Q 2 = 1.0 (GeV/c) 2 in 10 bins of 1.17 < W < 1.35 
GeV across the A resonance. A total of 16 independent response functions were extracted, of which 
12 were observed for the first time. Comparisons with recent model calculations show that re- 
sponse functions governed by real parts of interference products are determined relatively well near 
the physical mass, W = Ma ~ 1.232 GeV, but the variation among models is large for response 
functions governed by imaginary parts and for both increases rapidly with W > Ma- We per- 
formed a multipole analysis that adjusts suitable subsets of l-n < 2 amplitudes with higher partial 
waves constrained by baseline models. This analysis provides both real and imaginary parts. The 
fitted multipole amplitudes are nearly model-independent — there is very little sensitivity to the 
choice of baseline model or truncation scheme. By contrast, truncation errors in the traditional 
Legendre analysis of N — > A quadrupole ratios are not negligible. Parabolic fits to the W depen- 
dence around Ma for the multiple analysis gives values for Re(Si+/Mi+) = (—6.61 ± 0.18)% and 
Re(£a + /Mi+) = (-2.87 ± 0.19)% for the p -k° channel at W = 1.232 GeV and Q 2 = 1.0 (GeV/c) 2 
that are distinctly larger than those from the Legendre analysis of the same data. Similarly, the 
multipole analysis gives Re(5'o+/A/i+) = (+7.1 ± 0.8)% at W = 1.232 GeV, consistent with recent 
models, while the traditional Legendre analysis gives the opposite sign because its truncation errors 
are quite severe. Finally, using a unitary isobar model (UIM), we find that excitation of the Roper 
resonance is dominantly longitudinal with p S 1/2 = (0.05 ± 0.01) GeV" 1/2 at Q 2 = 1.0 (GeV/c) 2 . 
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The ReS'o+ and Re_Bo+ multipoles favor pseudovector irNN coupling over pseudoscalar coupling or 
a recently proposed mixed-coupling scheme, but the UIM does not reproduce the imaginary parts 
of 0+ multipoles, especially Im£o+, well. 

PACS numbers: 14.20.Gk,13.60.Le,13.40.Gp,13.88.+e 

I. INTRODUCTION 

The electroexcitation of a nucleon resonance is described by two transverse form factors, Ai/ 2 {Q 2 ) and A 3 / 2 (Q 2 ), 
and one longitudinal or scalar form factor, S±/2(Q ), where the subscript denotes the projection of the resonance spin 
upon the virtual photon momentum and the dependence upon photon virtuality, Q 2 , describes the spatial structure 
of the transition. Measurement of all three transition form factors could provide stimulating tests of QCD-inspired 
models of baryon structure For example, one might be able to determine the admixtures of small components 
with mixed symmetry or orbital excitation into wave functions dominated by SU(6) symmetry. Alternatively, it has 
been speculated that the Roper resonance at 1440 MeV could be a hybrid baryon based upon gluonic excitation 
0. The electromagnetic transition form factors would be distinctly different for a hybrid baryon or for a radial 
single-quark excitation However, the electroexcitation of a resonance, R, by a virtual photon, j v , in a reaction of 
the form j v N — > R — > Nx, where x could be one or even several mesons, inevitably is accompanied by nonresonant 
contributions and is complicated by final-state interactions. Furthermore, the nucleon resonances are broad and 
overlapping. Therefore, measurements of differential cross sections alone are not sufficient to permit clean, model- 
independent determination of the multipole amplitudes for meson electroproduction or, ultimately, the transition form 
factors. 

There has been long-standing interest in deformed components of the N and A wave functions The dominant 
amplitude for pion electroproduction at the A resonance is the M\ + amplitude, but smaller S\+ and E\ + amplitudes 
arise from configuration mixing within the quark core [f| , often described as quadrupole deformation, or from meson 
and gluon exchange currents between quarks 6], or coupling to the pion cloud outside the quark core 0,0- Thus, 

(3/2) 

the quadrupole deformation is related to the ratios of isospin-3/2 electroproduction multipole amplitudes R EM = 

Re(^+ /2) /Af^ /2) ) and R { g[ 2) = Re(5^ /2) /M 1 ( + /2) ) evaluated at W = M A where M A = 1.232 GeV is the physical 
mass. These quantities are often labelled as EMR and SMR instead and we will use both notations interchangeably. 
It has been argued that the intrinsic N and A quadrupole moments arc small 9] and that the observed EMR and 
SMR are dominated by nonvalence degrees of freedom. 

Most previous determinations of EMR and SMR analyzed differential cross section data using Legendre expansions 
truncated according to the assumption of M\ + dominance. The data prior to 1990 generally suffered from poor 
statistics, limited acceptance, and relatively large systematic uncertainties, leaving even the sign of EMR at low Q 2 
undetermined ^(j. More recent experiments using polarized photons pH Il2| now find Rem ~ —2.5% at Q 2 = 
with relatively little model dependence p|. Similarly, recoil polarization for parallel kinematics gave Rsm ~ —6.4% 
for Q 2 = 0.12 (GeV/c) 2 0. Furthermore, using more pre cise measurements of the azimuthal dependence of the 
differential cross section give Rsm ~ —6.1% at Q 2 — > |l5l lili. Il7j and have mapped the Q 2 dependence up to 4 
(GeV/c) 2 [ll[l!II3, but these analyses still rely upon Mi + dominance. However, there are several indications that 
this truncation may be inadequate. First, at Q 2 » 0.13 (GeV/c) 2 there is a strong disagreement between the the SMR 
value obtained at Bonn [2l| detecting a forward pion and those obtained at MIT-Bates and Mainz [T^J detecting 
a forward proton that might be explained by an unexpectedly large 5o+ amplitude |22(. Second, the dynamical 
models fail to reproduce the induced polarization for parallel kinematics at Q 2 < 0.2 (GeV/c) 2 where the pion cloud 
is important [23, 0|- Similar problems have also been observed at Q 2 = 0.4 and 0.65 (GeV/c) 2 [25|. Third, it 
appears to be difficult to obtain a consistent description of the real and imaginary parts of the longitudinal-transverse 
interference for Q 2 < 0.25 (GeV/c) 2 and W < 1.23 GeV where one might expect M 1+ dominance to suffice |l6l |26| . 
but a multipole analysis is not possible without more complete angular coverage and sensitivity to phases. 

In principle, this model dependence can be reduced considerably by measurement of recoil and/or target pol arization 
observables that are sensitive to the relative phase between resonant and nonresonant mechanisms |27ll28ll29ll30| . Po- 
larization often enhances the sensitivity to small amplitudes by interference with a large amplitude. Helicity-dependent 
recoil polarization is sensitive to real parts and helicity-independent recoil polarization is sensitive to imaginary parts 
of products of multipole amplitudes. Furthermore, by measuring the azimuthal dependencies one can determine re- 
sponse functions sensitive to the linear polarization of the virtual photon. Thus, it may be possible to perform a 
nearly model-independent multipole analysis if complete angular distributions are measured for a comprehensive set 
of recoil-polarization observables. These multipole amplitudes can then be interpreted with the aid of dynamical or 
coupled-channels models that enforce unitary and relate the complex multipole amplitude to the real transition form 
factor. 
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Note that coarsely binned measurements of asymmetries with respect to longitudinal target polarization have been 
made recently by Biselli et al. |3l| for the jot channel and by De Vita et al. [32] for the mr + channel in the resonance 
region with 0.35 < Q 2 < 1.5 (GeV/c) 2 . However, those data have not been resolved into response functions and are 
not suitable for detailed multipole analysis. 

In this paper we report the first extensive measurements of angular distributions for recoil polarization in pion 
electroproduction. We have measured recoil polarization response functions for the p(e, e'p)ir° reaction at Q 2 s» 1 
(GeV/c) 2 near the A resonance, obtaining angular distributions for a total of 16 independent response functions in 10 
steps of W covering 1.17 < W < 1.35 GeV; the angular coverage and statistical precision are best in the central region, 
1.21 < W < 1.29 GeV. The data for W — 1.23 GeV were reported in Ref. |3^. Twelve of these response functions are 
observed here for the first time. We compare a traditional truncated Legendre analysis with a more general multipole 
analysis of these data. Although the two analyses are qualitatively consistent, the multipole analysis shows that 
truncation errors in the traditional Legendre analysis of the quadrupole ratios are not negligible. Furthermore, we 
find that excitation of the Roper resonance is dominantly longitudinal. 

Section ITT1 defines the response functions, Sec. IIIII describes the experiment, and Sections IIVI and Ivl describe the 
cross section and polarization analyses, respectively. We compare the results with selected models in Sec. IVI Al and 
present Legendre and multipole analyses in Sections IVIBI and IVI CI Further discussion of the relationship between 
the two analyses is given in Sec. IVIII An interpretation of the multipoles amplitudes using a unitary isobar model is 
then given in Sec. IVIIII Finally, Sec. IIXI summarizes our conclusions. 



II. RESPONSE FUNCTIONS 



A. Definitions 



The differential cross section for recoil polarization in the pion electroproduction reaction p(e, e'N)ir can be expressed 
in the form 

r 7 a [1 + hA + U-S] (1) 



dkrdVLpdVL 



where tr is the unpolarized cross section, h is the electron helicity, S is the spin direction for the recoil nucleon, A is 
the beam analyzing power, U = P + hP' is the recoil polarization, 

1 2ir 2 hQ 2 \-e 1 ' 



is the virtual photon flux for initial (final) electron momenta fej [kf), e = ^1 + 2|p tan 2 4fJ is the transverse 

polarization of the virtual photon, 9 e is the electron scattering angle, and fc 7 = (W 2 — m 2 )/2m p is the laboratory 
energy a real photon would need to excite the same transition. Note that electron kinematics and solid angle d£l e are 
normally expressed in the lab frame while hadron kinematics and nucleon solid angle dfl* are expressed in the irN 
cm frame. Figure ^ illustrates the kinematics of this reaction and the definitions for polarization vectors. 

It is convenient to express polarization vectors in a basis where £ is along the nucleon momentum in the cm frame, 
n oc q A i is normal to the reaction plane, and t = h A I. The azimuthal dependence can then be extracted and 
the observables can be decomposed into kinematical factors, v a , which depend only upon electron kinematics and 
response functions, R a {x, W, Q 2 ), which carry the hadronic information, such that 



a = u \vlRl + vtRt + vltRlt sin 6 cos <\> + vttRtt sin 2 cos 2(j)\ (3a) 

Aa = vq [z^ T i?^ T sin # sin </>] (3b) 

P t a = Vo [vltR-lt sin cf) + vttRtt sm $ sin 2(f\ (3c) 

P n a = u \(y L R n L + u t Rt) sin 9 + v LT R n LT cos <j> + v T tRtt sin cos 2<p] (3d) 

Pga — vq \vltR 1 lt sm ^ sm 4> + vttRtt sm2 ^ sm ^ < f > ] (^ e ) 

P[a = vq [v' lt R'lt cos 4> + v' TT R!^ T sin 9] (3f ) 

P' n a = vq [v'ltR'lt sin <P] (3g) 

P' t a = VQW LT R'l T shx9cos^ + v TT R^ T \ (3h) 
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where the response functions depend upon W, Q 2 , and x — cos 9 where 9 is the pion angle relative to q in the irN 
center of mass frame. The azimuthal angle <p also refers to the pion momentum. The overall factor v$ permits phase- 
space factors to be extracted from the response functions, thereby simplifying multipolc expansions. This factor is 
usually taken as the ratio between the cm momentum in the final state and the cm momentum a real photon needs 
for the same transition, such that vq = k/qo where 



2\2 



(W 2 + m%- m 2 ) 

m 2 



K (4a) 



W 2 -ml 

q = ^ 4b 

Regrettably, no single convention for the signs and normalizations of the response functions has gained wide accep- 
tance. We chose 

vt = 1 (5a) 

vtt = e (5b) 

v' TT = Vl-e 2 (5c) 

v L = e s (5d) 

v LT = ^/2e s {l + e) (5e) 



= ^2e s {l - e) (5f) 



where 



Q 2 rfi , 

6 S = —ft (6) 

q- 

and use the azimuthal angle for the pion in Eq. ©. Note that although e is invariant, es is not and that Eq. is 

evaluated in the ttN cm frame. When there is insufficient information to perform a Rosenbluth separation, we employ 
e-dependent combinations 

v t Rl+t = v L R L + v T R T (7a) 

vtRI+t = vlRI + vtRt (7b) 
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TABLE I: Truncated Legendre expansion of the differential cross section. Legendre coefficients based upon the assumption 
of Mi + dominance are obtained as either the real or imaginary parts of M* + B v , m . The final column lists the labels used in 
figures. 
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B. Legendre Expansions 

Legendre expansions often provide the most efficient representation of the angular dependence of a response function. 
Each of the response functions can be represented by a Legendre expansion of the form 

00 

R v (x, W, Q 2 ) = Q 2 )P m (x) (8) 

771 = 

where r\ represents all of the labels needed to identify a particular response function. Notice that by extracting 
the leading dependencies upon sin#, the response functions defined by Eq. © reduce to polynomials in x that are 
expected to be of low order for modest W and Q 2 . This convention should also improve the accuracy of acceptance 
averaging within bins of x. The Legendre coefficients A^W, Q 2 ) can be obtained by fitting response- function angular 
distributions for each (W,Q 2 ). Alternatively, often the most efficient method for extracting the x dependence of 
response functions for a particular (W, Q 2 ) bin is to perform a two-dimensional fit of the (x,(j>) dependencies of the 
appropriate observable, cross section or polarization, with the aid of the Legendre expansion. This type of analysis can 
be used to extrapolate response functions to parallel or antiparallel kinematics where interesting symmetry relations 
have been developed but where the experimental acceptance vanishes |34L l35| . 

Each of the Legendre coefficients can be further expanded in terms of products of pairs of multipole amplitudes, 
but these expansions quickly become unwieldy as the number of participating partial waves increases. Analyses of 
previous data for differential cross section have usually found that fits using truncated Legendre expansions based 
upon Mi + dominance describe the data well in the immediate vicinity of the A resonance. Tables lllllll display 
truncated expansions that are limited to s and p waves and that require Mi 4. to appear in every term. To streamline 
the notation, each coefficient is obtained from either the real or imaginary part of Mi + B vm where the phase type is 
indicated in the second column and B^ m is listed in the fourth column of these tables. The considerable redundancy 
among the B v ^ m functions can be used to perform a model-independent evaluation of the accuracy of Mi 4. dominance. 
The last column gives the response-function labels that will appear in multipanel figures. 

III. EXPERIMENT 

The experiment was performed in Hall A of the Continuous Electron Beam Accelerator Facility (CEBAF) at 
Jefferson Laboratory in three periods between May and August 2000 using standard equipment described in detail 
111 Ref. HI. Details of the design and performance of the focal-plane polarimeter (FPP) can be found in Ref. 37]. 
Further information specific to this experiment may be found in Refs. j3^,|3^,^3|- Here we summarize the salient 
features. 

The beam energy of 4531 ± 1 MeV was measured by two independent methods and the results were in agreement. 
The arc method uses wire scanners to measure the beam deflection angle between the accelerator and the experimental 
beam lines. The eP method measures electron and proton angles for elastic scattering from hydrogen using silicon 
strip detectors placed symmetrically in the vertical plane of a dedicated instrument. 

The beam current was monitored by a pair of resonant RF cavities that were periodically calibrated with respect to 
an Unser monitor. The Unser monitor is a Parametric Current Transformer that provides absolute current measure- 
ments and is calibrated with respect to a precision current source, but it suffers from a variable offset. The RF cavities 
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TABLE II: Truncated Legendre expansion of the helicity-dependent recoil polarization. Legendre coefficients based upon the 
assumption of M 1+ dominance are obtained as either the real or imaginary parts of M* + B Vtm . 
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TABLE III: Truncated Legendre expansion of the helicity-independent recoil polarization. Legendre coefficients based upon 
the assumption of Mi+ dominance are obtained as either the real or imaginary parts of Mi + B Virn . 
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are much more stable and can be used as continuous monitors. The calibration procedure alternates between beam 
off and current ramping between about 10 and 100 /iA in about 5 steps, using the beam off periods to determine the 
drift of the Unser offset and the current steps to determine the RF gains. This procedure was repeated several times 
durin g th e experiment, with negligible differences between calibrations. The accuracy of these monitors is better than 
0.5% |3|. 

The beam position and direction were measured by Beam Position Monitors (BPMs) located 7.524 and 1.286 m 
upstream of the target. Each BPM is a 4- wire antenna array tuned to the fundamental RF frequency of the beam 
and provides the relative position to within 100 /im. The absolute positions of the BPMs are calibrated with respect 
to wire scanners that are surveyed regularly. 

The polarized electron beam was produced by circularly polarized light from a 780 nm diode laser, pulsed at 1497 
MHz, illuminating a strained gallium arsenide wafer. A Pockels cell was used to reverse the laser polarization at 30 Hz. 
The beam polarization was measured nearly continuously using a Compton polarimeter, with systematic uncertainties 
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FIG. 2: (Color online) Beam polarization measurements: gray boxes indicate Moller measurements and points Compton 
measurements, with systematic uncertainties indicate for both. Vertical dashed lines mark movements of the laser spot. The 
red vertical line makes the change of polarized source. 



estimated to be about 1% |40j |. In addition, periodic measurements were also made with a Moller polarimeter, with 
systematic uncertainties of about 2.4%. Statistical uncertainties were negligible for both methods. Fig. |21 summarizes 
the polarization measurements, where the vertical dashed lines indicate movement of the laser spot needed to maintain 
high current. The polarization averaged about 72% for the first two running periods, but for the third it was necessary 
to use a different source with lower polarization, about 65%. Fig. shows that the time dependence of the beam 
polarization is minimal between spot changes. Moller results were used during the few occasions when the Compton 
measurements were unavailable. A comparison between the five beam polarimeters at Jefferson Lab 41] indicates 
that Compton measurements in Hall A are consistent with Mott measurements at the injector and that the ratio 
between Hall A Moller and Compton analyzing powers is approximately 1.034 ± 0.028, which is consistent with the 
observation in Fig.[21that the Moller results are generally slightly larger than the Compton results. However, because 
the statistical precision of this ratio remains poor and relatively few runs rely on the Moller polarimeter, a correction 
for the relative normalization of the two polarimeters was not applied. 

The production target was a cryogenic loop for liquid hydrogen operating at 19 K and 0.17 MPa with a density of 
0.723 g/cm 3 . The beam passes through a cylindrical aluminum cell that is 15 cm long and 6.35 cm in diameter with 
sidewall thickness of 178 /im and entrance and exit windows approximately 71 and 102 fim thick, respectively. To 
avoid damage to the target cell and to minimize luminosity variations due to local heating or boiling, the beam was 
rastered over a 4 x 4 mm 2 pattern. Measurements of rates versus current demonstrate that the average target density 
decreases by less than 2% for rastered currents in the range 20-60 /iA relevant to the cross section measurements [3!j . 
Furthermore, the electron singles rate provided a continuous luminosity monitor because the electron spectrometer 
remained fixed throughout the experiment. 

Additional data were taken with about 5 ^A of unrastered beam on a thin carbon foil and on dummy targets 
consisting of aluminum foils separated by either 4, 10, or 15 cm that mimic the windows of a cryogenic cell. These 
data were used to determine the mispointing of the spectrometers and the data for the 15 cm dummy cell were used 
to determine the contribution of the cell walls for the production data. 

Particles were detected by a pair of high resolution QQDQ spectrometers of identical design, with electrons observed 
in the left and protons in the right spectrometer. These spectrometers have a central bend angle of 45° and nominal 
acceptances of ±4.5% in momentum, ±60 mr in vertical angle, and ±30 mr in horizontal angle. The cross section data 
taken during the first of three running periods were taken with 80 mm tungsten collimators at distances of 1.109 m 
for the left or 1.100 m for the right with nominal apertures of 6 msr. Polarization data taken in the next two running 
periods were acquired with open apertures. The entrance windows were 178 /im Kapton while the exit windows were 
100 fim titanium. 

The position of the electron spectrometer was determined by survey. The initial position of the proton spectrometer 
was also determined by survey for each run period and subsequent angles were measured using accurately placed floor 
marks observed through a closed-circuit television camera mounted on a linear translation stage without parallax. 



8 



These raw angle measurements were then corrected for roll and pitch as measured by bi-axial inclinometers. However, 
because the spectrometers do not rotate about a fixed pivot, it is necessary to correct for possible mispointing. In 
principle, mispointing can be determined by a system that uses a linear variable differential transformer (LVDT) 
to measure the gap between an arm parallel to the spectrometer mid-plane and the outer surface of the scattering 
chamber, but this method proved to be unreliable. The LVDT mispointing data were compared with mispointing 
measurements deduced from scattering data for the carbon target. The latter were found to be more stable and were 
used in the subsequent analysis. The spectrometer offsets deduced for the electron arm, which was stationary for each 
run period, were consistent with constant values that were less than 2 mm and consistent with survey. The offsets for 
the hadron arm varied between about -4 and +1.5 mm with the motion of the spectrometer, but were reproducible 
without motion [38l| . 

Both arms use two vertical drift chambers (VDC) for tracking that are inclined by ±45° with respect to the central 
trajectory and are separated by 50 cm. Each chamber contains two planes of sense wires inclined at ±45° with 
respect to the midplane (uv configuration). Valid tracks typically produce signals on 3-5 sense wires. These detectors 
typically provide position resolutions of a x ,y = 100 /im and angular resolutions of crg^ — 0.5 mrad. Further details 
about the VDCs are provided in Ref. |4^ . 

Both arms also use two trigger planes, SI and S2, consisting of 6 plastic scintillator paddles that are 5 mm thick, 
overlap by about 0.5 cm, and are viewed by photomultiplier tubes on both ends. Five types of trigger were produced 
by the trigger supervisor module using signals from the two trigger planes in each spectrometer. A Tl (T3) trigger 
in the electron (proton) arm requires coincidence between paddle i of SI and paddle j of S2 with \i — j\ < 1, with 
PMT signals above threshold in both ends of both paddles. A T2 (T4) trigger in the electron (proton) arm misses 
one or more of the PMT signals or fails the directivity requirement. The T5 primary coincidence trigger requires Tl 
and T3 within 100 ns of each other, where T3 and T4 include adjustable delays depending upon hadron momentum. 
The time resolution for the coincidence trigger is typically about a = 0.3 ns. Events corresponding to any of these 
triggers can be recorded, but triggers other than T5 were generally prescaled to limit computer deadtime. 

The electron arm also contained a CO2 gas Cerenkov and lead-glass pre-shower and shower detectors for particle 
identification, but these were not used for this experiment because electrons were overwhelmingly dominant for these 
conditions. A simple cut on flight time between scintillator planes was sufficient. The gas Cerenkov detector normally 
present in the proton arm was removed for this experiment. For the first running period a single large scintillator 
paddle, labeled SO, was included in the proton trigger. After the first running period the SO scintillator was replaced 
with an aerogel Cerenkov detector with n = 1.025 and the aerogel signal was used in the proton trigger to suppress pion 
background with momentum greater than 0.63 GeV/c. We could then make polarization measurements at forward 
angles with higher beam currents without excessive deadtime. 

Recoil polarization measurements were made using the focal plane polarimeter (FPP) installed in the proton arm. 
More complete details of the construction and operation of the FPP can be found in Refs. [36ll37| . The FPP consists 
of two front straw chambers, a carbon analyzer of adjustable thickness, and two rear straw chambers. The front 
chambers were separated by 114 cm and each contains 6 planes in a vvvuuu configuration. The rear chambers were 
separated by 38 cm. Chamber 3 has a uuvvxx and chamber 4 a uuuvvv configuration, where x denotes sensitivity to 
the dispersive direction. The analyzer consisted of five graphite plates with thicknesses of 1.9, 3.8, 7.6, 15.2, and 22.9 
cm separated by about 1.6 cm that can be deployed in any combination to optimize the analyzing power for specified 
proton momentum. 

Spectrometer settings were chosen to obtain (W,Q 2 ) — (1.232 GeV, 1.0 (GeV/c) 2 ) for a nominal beam energy of 
4.535 GeV. Thus, the electron spectrometer remained fixed at 14.095° with a central momentum of 3.660 GeV/c. The 
nominal settings for the proton spectrometer, with accumulated charge, are summarized in Table llVl It is convenient 
to define 8 pq as the nominal center-of-mass angle between the nucleon momentum and the momentum transfer in 
the lab frame with positive (negative) sign corresponding to forward (backward) angles with respect to the beam 
direction. The experiment was divided into three periods. The first period was used to scan the angular distribution 
and to obtain differential cross sections for all settings except 9 pq = 180° using the 6 msr collimator in the hadron 
arm and relatively low currents to limit deadtime. The second and third periods used the highest possible currents to 
obtain adequate statistics for recoil polarization and used the aerogel detector for forward angles to suppress deadtime 
due to accidental tt + coincidences. 

Electrons were identified using time of flight between the SI and S2 scintillator planes. Protons were identified using 
the velocity measured by time of flight between the SI and S2 scintillator planes and the geometric mean between 
pulse heights for the two sides of SI. Fig. [3] shows a two-dimensional distribution for 9 pq = —155° where the pion 
background is relatively large. The polygonal cut selects protons with little loss and with little contamination. A cut 
was also made on the reconstructed vertex position. This cut was typically ±2 cm; the resolution is limited mostly 
by the forward angle of the electron spectrometer. Finally, a cut on the missing energy and momentum, (E m ,p m ), 
for the pir° final state was used in the polarization analysis to suppress background from elastic scattering. This 
background was visible only for the pq = —50° and —90° settings; the pw event selection is illustrated in Fig. 0]for 
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TABLE IV: Settings for the proton spectrometer. Here 9 P is the laboratory angle with respect to the beam, p p is the central 
momentum, and 6 pq is the cm angle with respect to q. The fifth column lists the 12 C thickness used for the FPP and the final 
column lists the nominal precession angle for the central momentum. 
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FIG. 3: (Color online) Proton identification using the correlation between velocity and energy deposition in SI. The polygon 
selects protons. Notice that the intensity scale is logarithmic. 



9 pq — —50°. However, this cut was not needed for the cross-section analysis because the more restrictive acceptance 
cuts described in Sec. IIV Al eliminated the elastic background quite effectively. 
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FIG. 4: (Color online) Pion production events for 6 vq = —50° are selected by the polygon in (Em,Pm)- The diagonal population 
arises from the elastic radiative tail. The intensity scale is logarithmic. 



IV. CROSS SECTION ANALYSIS 



A. Acceptance 

In each spectrometer a track is defined by the four reconstructed target variables (S, y, </>)tgt! these variables are 
defined according to TRANSPORT conventions. Rather than attempt to visualize the four-dimensional volumes 
populated by event coordinates for each spectrometer, one normally inspects two-dimensional projections, such as 
those shown in Fig. [5] for the electron spectrometer. An event that is safely near the center of the distribution in 
one projection may be found near the edge of another where either the experimental acceptance or the calibration 
of the magnetic optics may be suboptimal. Therefore, it is useful to construct a measure of the distance between 
event coordinates and the boundary of a multidimensional acceptance volume that is based upon two-dimensional 
projections that are more amenable to visualization. We employ a variation of the i?-function method that was 
originally developed by V. Rvachev 0, and applied to (e, e'p) reactions by M. Rvachev pjj. |47| . We begin 
by drawing regular polygons that encompass most of the pion electroproduction events within each of the 6 two- 
dimensional projections formed by pairs of variables. The two coordinates (x, y) are then rescaled according to 



where the denominators represent the extremes of the polygonal boundaries and where TRANSPORT conventions ensure 
that the polygon is centered near the origin. For each event with coordinates (x',y') in projection i, we define di to 
be plus or minus the minimum distance to a side of the polygon for that projection, with a positive (negative) sign 
used for points that are inside (outside) the enclosed area. Figure illustrates these signed distances. Differences 
between the units employed for the various transport variables are then compensated by defining a normalized distance 
& = di/^/Al where A{ is the area enclosed by polygon i in terms of normalized coordinates. Thus, & represents the 
normalized distance of an event from the acceptance boundary in projection i with positive signs inside, negative 
signs outside, and zero on the boundary. The inner region of Fig. 0with £j > 0.05 represents events that are safely 
inside the polygonal boundary represented by £j = 0. Finally, for each spectrometer we define R = ±min{£j} with 
a positive sign when all £j are positive and a negative sign whenever an event falls outside the border for any of the 
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FIG. 5: (Color online) Polygonal boundaries for two pairs of reconstructed electron variables for coincidence events. 
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FIG. 6: Definition of signed distances from the boundary of a two-dimensional projection. 



projections. 

This method can be applied to each spectrometer independently or to coincidence data for which kinematic cor- 
relations between tracks in the two spectrometers suppresses accidentals. For the coincidence method we define 
R = min(i? e , R p ) as the smaller of the R- functions for the two spectrometers. However, the disadvantage of the co- 
incidence method is that the polygons must be constructed for each kinematical setting independently. Nevertheless, 
we chose the coincidence method for this experiment because there were often substantial accidental populations. 
Further details concerning the algorithms for the present implementation of the i?-function method can be found in 
Ref. lal. 



B. Simulation 



The differential cross sections for each kinematical bin were obtained from comparison between experimental and 
simulated yields. The simulations used MCEEP ji^l to apply radiative corrections to a theoretical model, resulting in 
a six-fold differential cross section, and to integrate the theoretical cross section over the experimental acceptance. 
MCEEP samples the phase space uniformly, over a volume larger than the experimental acceptance, and evaluates the 
yield 

Yi = C J K,i<fa <g> K (10) 

where C is the luminosity, JCi represents the acceptance function for bin i, d 6 <r represents the model cross section for 
each event, ® represents convolution, and 1Z represents resolution functions for quantities measured in the focal-plane. 
Here /Q = 1 if R > R C ut or otherwise based upon target variables that are reconstructed from focal-plane coordinates 
convoluted with resolution functions. 
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FIG. 7: Regions selected by an J?-cut in two dimensions. 



The model cross section was based upon tabulated multipole amplitudes for MAID2000 The kinematics for 
each event were used to interpolate the multipole amplitudes with respect to (W, Q 2 ) and then to compute the 5-fold 
differential cross section in the laboratory frame. 

The radiative corrections include bremsstrahlung in the target before and after scattering, internal soft-photon 
processes according to the Schwinger prescription, and radiation of hard photons using the Borie-Drechsel |50j pre- 
scription with the peaking approximation. Multiple scattering within the target and windows is included also. These 
corrections do not account for polarization effects. Further details can be found in Refs. [39l l48|| . 

Figure |H1 compares simulated and experimental distributions for target variables at 9 cm = using an acceptance 
cut R > 0.05 for both, and applying a common normalization factor to the simulation. Although a slight discrepancy 
can be observed for 9 tg in the electron spectrometer, the simulation reproduces the experimental distributions very 
well. Similarly, Fig. shows that the simulation also reproduces the distributions for reaction kinematics very well. 
Note that the experimental acceptance function, R e h, is shown without applying the cut. The optimal choice for i? cu t 
is somewhere below the center of the plateau in the ratio between experimental and simulated yields. The systematic 
uncertainty due to the acceptance function, estimated from the flatness of the plateau, is typically less than 1%. 

The missing mass for this reaction is quite sensitive to laboratory angles; for example, at 9 pq = —155° the sensitivity 
to electron angle is 13 MeV/degree for our kinematics. Thus, comparing the simulation with data for 8 pq — —155° 
we adjusted the electron angle by —0.05°, which is well within the survey uncertainty, and then find good agreement 
with the missing mass peaks for all other settings as well. Furthermore, the width of the missing mass peak is 
underestimated by the simulation unless resolution functions are applied to the track coordinates. We found that 
Gaussian resolutions with a — 0.5 mm applied to the hit positions in each VDC plane provide good agreement with 
those widths. 



C. Background 

The acceptance cut R > 0.05 suppresses the background from the elastic radiative tail quite strongly. The residual 
contribution is less than 0.4% at 9 pq — —90° and much smaller at 9 pq = —50°. Hence, no corrections for this 
background were needed for the cross section or the unpolarized response functions. The background due to accidental 
coincidences was subtracted using time windows on both sides of the coincidence peak, applying the same particle 
identification and acceptance tests, and normalizing by width. 

D. Cross section calculation 

The virtual photoproduction cross section a for a particular kinematical bin was determined by scaling the model 
cross section <7 mo dei for that bin, evaluated for bin-centered kinematics, and applying various deadtime and efficiency 
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FIG. 8: (Color online) Comparison between observed (solid blue) and simulated (dotted red) distributions of target variables 
for # pq = with i? C ut > 0.05. E-arm refers to the electron and H-arm to the hadron spectrometer. 



corrections according to 

Y _ /cDT/EDT/abs M 1 « 

<7 = T? Cmodcl (11) 

'MC EtriggcrCtrack 

where Y and Ymc are the observed and simulated yields, /cdt corrects for computer deadtime, /edt corrects for 
electronics deadtime, cogger corrects for trigger efficiency, etrack corrects for wire chamber and tracking efficiency, and 
/abs corrects for proton absorption in materials between the scattering and detection. No correction was made for 
variation of luminosity with current because no systematic variation was observed in the luminosity monitor for the 
currents employed in the cross section measurements (20-60 /LtA). 

The computer deadtime was determined by comparing the coincidence scaler with the number of coincidence 
events recorded. In addition, the trigger supervisor has an internal deadtime of approximately 100 ns, such that the 
electronics deadtime for a 1 MHz rate in the trigger scintillators is about 10%. The electronics deadtime was measured 
by sending pulser signals to one scintillator paddle in each arm and comparing the number of pulser signals recorded 
with the number counted by a scaler, correcting for computer deadtime. The dependence of the electronic deadtime 
upon strobe rate was then parametrized. The systematic uncertainty in the correction for electronic deadtime was 
estimated to be about 1% at the highest rates |51j. 

The event reconstruction software rejects events with more than one track in either spectrometer. For the electron 
spectrometer 10-12% of the events contained multiple tracks while for the proton spectrometer 1-12% contained 
multiple tracks depending upon the momentum and angle settings. We assume that the fraction of multiple-track 
events that contain a particle that would have satisfied the particle-identification criteria and other tests is the same as 
that for single-track events and apply corrections for each arm independently. In addition, we required valid tracks to 
contain 3-8 hits in each VDC plane. For the two settings with significant population by elastic scattering, the elastic 
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FIG. 9: (Color online) Comparison between observed (solid blue) and simulated (dotted red) distributions of kinematical 
variables for 6> pq = with -R cu t > 0.05. The lower left panel shows the uncut distribution of R e h- 



scattering events were excluded from the calculation of tracking efficiency to minimize position-dependent effects upon 
trigger efficiency and to improve factorization of the tracking efficiencies for the two arms. 

The triggers in each arm require coincidence between two scintillator planes and test the track direction. Thus, 
the trigger efficiency compares the total number of valid triggers with the total number of events with a least one 
hit in a scintillator. For the electron arm, we require events in both the numerator and the denominator to satisfy 
the Cerenkov test for electrons, and to contain only one track. For the proton arm we also use one-track events but 
use the SO scintillator instead of the Cerenkov detector. The net trigger efficiency of approximately 96.7%, with a 
systematic uncertainty of about 1%, is then the product of the efficiencies for the two arms 

Finally, we used a compilation of proton reaction cross sections to estimate the probability for proton absorp- 
tion between scattering and detection. The net correction factor varied between 1.008 and 1.017 depending upon 
momentum. 

E. Cross section results 

We assume that the ratio between observed and simulated yields over the acceptance for a kinematic bin is very 
nearly equal to the ratio between actual and model differential cross sections for the central kinematics of that bin. 
The accuracy of this assumption depends, of course, upon the bin size and the curvature of the differential cross 
section with respect to the binned variables. Events were accepted for Q 2 = 1.0 ± 0.2 (GeV/c) 2 . We used 10 bins in 
W between 1.17 and 1.35 GeV of width ±0.01 GeV, 20 bins of x between -0.95 and +0.95 in steps of 0.1, and 12 bins 
of <fi of 30° width. After dropping bins with negligible acceptance, approximately 1140 data were obtained for both 
differential cross section and beam analyzing power. These data are reported for central kinematics. 
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FIG. 10: Fits to the <j> dependence of unpolarized cross section data with (W, Q 2 ) = (1.23, 1.0). Each panel is labeled by the 
central x. Solid curves fit Legendre coefficients to the entire data set shown while dashed curves fit response functions within 
each panel independently. 

Figure ITU1 shows the dependence of the differential cross section for each x bin with (W,Q 2 ) = (1.23,1.0). 
The dashed curves fit Rl+t, Rlt, and Rtt for each (x, W, Q 2 ) independently to the dependence of Eq. l|5a|l . 
Unfortunately, this procedure did not permit model-independent separation of Rtt from Rl+t for x ~ because 
correlations were too large given the present acceptance. The solid curves fit Legendre coefficients to the (x, 0) 
dependence, thereby imposing a smooth x dependence that is not required by the extraction of unpolarized response 
functions. Nevertheless, both methods fit the data well and agree within the uncertainties estimated from covariances. 
Similar figures can be made for each (W 7 Q 2 ) bin, but are too numerous to display here. 

The Legendre coefficients fit to the unpolarized cross section for Q 2 = 1.0 (GeV/c) 2 are compared in Fig. 1111 with 
expansion coefficients for calculations based upon the MAID2003, DMT, SAID, and SL models obtained by inversion 
of Eq. (jSJ). (More details about model calculations are given in Sec. IVI A|) . Although these calculations suggest that 
the sp truncation is probably adequate in the immediate vicinity of the A resonance, it appears that additional terms 
may be necessary elsewhere. Therefore, in addition to fits based upon sp truncation, we show fits with one additional 
free parameter for each response function within the the central W range where the angular coverage and statistical 
precision are best. The models reproduce the even L + T coefficients relatively well, although the W dependence of 
the SAID calculation for A^ +T is somewhat too flat. The models also reproduce the low-order coefficients for Rlt 
and Rtt relatively well. For Rtt the additional coefficient is determined relatively well near the middle of the W 
range and is consistent in both sign and magnitude with most model calculations. The resulting curvature in Rtt 
is small but definitely visible. Similarly, the data are consistent with the small negative linear coefficient predicted 
for Rtt but cannot determine higher-order coefficients. The additional term for Rl+t appears to be rather weak. 
The extra terms have very small effects upon the fitted value for lower coefficients of the same parity, but negligible 
effect upon those of opposite parity. Note that A^ +T is appreciably stronger than MAID, DMT, or SL predictions 
and exhibits an upturn for W > 1.3 GeV that is absent from those models and that this result is not affected by the 
inclusion of terms beyond sp truncation. 

Figure ITTI also shows similar results obtained by Joo et al. [Lsl l52]| at Jefferson Laboratory using CLAS. Here we 
show their results for the higher beam energy, 2.445 GeV, that has better statistical precision. However, the two 
experiments used different binnings with respect to Q 2 . For the purposes of this comparison, we assume that the 
Legendre coefficients are proportional to the square of a dipole form factor and rescale the CLAS data for Q 2 = 0.9 
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FIG. 11: (Color online) Legendre coefficients fit to differential cross sections at Q 2 = 1.0 (GeV/c) 2 . These quantities are denned 
in Eq. JHJ where superscripts identify the response function and subscripts the degree of the Legendre polynomial. The filled 
triangles use the sp truncation while in the central W range the filled circles include an extra term for each response function. 
CLAS results scaled to Q 2 = 1.0 (GeV/c) 2 using a dipole form factor are shown as open green circles for Q 2 = 0.9 (GeV/c) 2 
and open green squares for Q 2 = 1.15 (GeV/c) 2 ; only data for a beam energy of 2.445 GeV are shown. These results are 
compared with MAID2003 (red solid), DMT (green dashed), SAID (blue dash-dotted), and SL (cyan dotted) calculations. 

and 1.15 (GeV/c) 2 to a common value of 1.0 (GeV/c) 2 . Note, however, that for a given W, e is higher for Q 2 = 0.9 
than for 1.15 (GeV/c) 2 and that e for our experiment is higher than for either of the CLAS data sets. We observe 
good qualitative agreement between these data sets, but there are significant differences in detail. For example, our 
A% +T is systematically stronger for low W than in CLAS data. Nor does the form-factor scaling prescription bring 
the two CLAS data sets for LT coefficients in agreement with each other, but the higher Q 2 data also appear to show 
more scatter. On the other hand, the curvature we see in the x dependence of Rlt clearly requires at least one term 
beyond sp truncation; this is shown in more detail in Sec. IVIBI Perhaps the omission of A\ T from the CLAS analysis 
is partly responsible for discrepancies in the lower coefficients. 

V. POLARIZATION ANALYSIS 
A. Polarization Analysis using Likelihood Method 

Let T = (T t ,T n ,Ti) = P + hP e P' represent the proton recoil polarization at the target in the irN center of mass 
system, where h denotes the sign of the electron helicity and P e is the magnitude of the beam polarization, and 
let F = (F x ,F y ,F z ) represent the polarization at the focal-plane polarimeter with z along the nucleon momentum 
and y leftward with respect to the vertical plane containing the nucleon momentum. These vectors are related by 
a spin transport matrix S, representing a sequence of transformations from the target cm frame to the local FPP 
coordinate system, such that F — ST. The spin transport matrix is evaluated for each event. Details of the individual 
transformations are given in Appendix lAl 

The polarization components at the target were extracted from the azimuthal distribution for scattering by the 
FPP analyzer using the method of maximum likelihood |53j |. The likelihood function takes the form 

L = II o~ ( 1 + £ -e * sin 0ft>p + e i/ cos 0ft>p) (12) 

events 

of a product of the scattering probabilities for each event that satisfies the selection criteria for a given kinematical 
bin. The azimuthal scattering angle </>f pp is measured counterclockwise from the final x axis for each event and £ 
represents the false (instrumental) asymmetry, discussed in Sec. IV El The e coefficients are given by 



e a — A y (9{ pp ) ^ SapTp 





(13) 
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where A y (8f pp ) is the analyzing power for polar scattering angle 0{ pp , a € {x,y, z} identifies the polarization compo- 
nents at the FPP, (3 e {t, n, £} identifies components of T at the target in the irN cm frame, and S a p are elements 
of the spin-transport matrix given by Eq. (|A1|) . Although the scattering probability for each event is independent of 
the longitudinal polarization, the variation of spin transport within the experimental acceptance offers access to all 
three components of polarization at the target. 

An iterative method for extracting the target polarization that maximize the likelihood is presented in Appendix 
IbI If the asymmetries (e x , £ y ,£) are small, the problem reduces to the linear system 

V = A-R (14) 

where 

R=(P t ,P n ,P t ,Pl,P^,PD (15) 

is the result vector, 

v -?iT5 (16) 

is an element of the measurement vector, and 

is an element of the design matrix where the Greek indices {a, j3 = 1, 6} correspond to elements of the result vector 
and the Latin index i enumerates events that satisfy the selection criteria for a particular kincmatical bin. Elements 
of the result vector represent acceptance-averaged components of recoil polarization that are taken to be constant 
within each kinematical bin. Conversely, the elements of the measurement vector and design matrix accumulate 
contributions 

£ = a sin (/> {pp + 6 cos <f> ipv + c sin 2</> fpp + d cos 2</> fpp 
Ai = A(9 lpp )(S y t cos^fpp - S xt sm(/>{p P ) 
X-2 = A(9f pp )(S yn cos 4>{ pp - S xn sin<f>f pp ) 
A 3 = A{9i pp ){Syi cos0f pp - S'rfsin^fpp) 
A4 = hP e X\ 
A5 = hP e \i 
As = hP e \s 

that are evaluated independently for each event, where the event indices have been suppressed. 

B. Extraction of polarized response functions using Likelihood Method 

Binning with respect to 4> can be avoided by using Eq. @ to express the likelihood function 

l = n ^( i +t+v-R) as) 

events 

in terms of response functions 

n ( Tjt pt nil nn nil jyl r>l jylt r>lt nln p/£ r>lt \ 

it — [n LT: n TT , n L+T , n LT , jtl tt , ix lt , tt TT , n LT , ri TT , n LT , n LT , n TT ) yij) 
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FIG. 12: (Color online) Kinematical acceptance and the binning for the polarization analysis. Left: angular acceptance for 
9n,4>n in the cm frame. Right: (W, Q 2 ) acceptance. 

with coefficients 

ar\\ = Xivqvlt sin</> 

a T)2 — XiVqVtt sin (9 sin2</> 

a- 774 = \2V0VLT cos <j) 

arj 5 — \2VqVtt sin cos 20 

ar] e — A3 i>q Vlt sin sin <p 

arj-j — X^vqVtt sin 2 9 sin 20 

(7 778 = A/t^o^r cos 
(T r/9 = \^vqv' tt sin 

5" ?7l0 = A5^o^LT s i n '/ , 

ct 7711 = A6 vq v' lt sin cos <f> 

5" 7712 = \VQV' TT 

that incorporate the azimuthal dependencies event by event. The coefficients for response functions depend upon 
the unpolarized differential cross section for each event, which varies within the kinematical bin. This cross section 
was obtained by scaling the model cross section (MAID2000) calculated at the event kinematics by the ratio between 
the sp Legendre fit to the experimental cross section and the model cross section for the central kinematics of the 
bin. This Legendre fit is discussed in Sec. IIV El The Legendre parametrization can sometimes produce nonpositive 
cross sections for some events with kinematics at the edges of the acceptance; for those events we simply use the 
MAID2000 cross section and recognize that these extreme kinematics contribute very little to acceptance-averaged 
quantities anyway. 

C. Track reconstruction and selection 

The chambers were aligned with respect to each other and the VDCs using straight-through events obtained by 
removing the carbon analyzer. The track reconstruction algorithms are described in Ref. 38]. For our purposes it is 
sufficient to note that the hit multiplicity within the straw chambers is sufficient to define tracks before and after the 
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carbon analyzer. Thus, we can impose a requirement that the scattering vertex lie within the carbon plates used for 
a particular measurement. We also require that the polar scattering angle be in the range 5° < 8f pp < 20°, where the 
lower limit enhances the analyzing power by suppressing unpolarized Coulomb scattering and the upper limit keeps 
instrumental asymmetries small. Finally, to minimize false asymmetries due to the finite size of the rear chambers, we 
impose a cone test that demands that the entire cone subtended by the polar scattering angle for each track intercepts 
both rear chambers. The rear chambers are actually large enough that only a few percent of the events in the accepted 
6*f pp range fail the cone test. 



D. Calibration 



We fitted an extension of the McNaughton parametrization [54j of the p+ 12 C analyzing power using earlier Hall 
A data supplemented by new measurements of elastic scattering by the proton for momenta of 0.818, 1.066, 1.188, 
and 1.378 GeV/c in order to provide analyzing power data closer to some of the present kinematical settings. Our 
measurements of Ge p /Gmp at Q 2 = 1-0 and 1.4 (GeV/c) 2 are in good agreement with those of Ref. [55| . 



E. False asymmetry 

The one-photon exchange approximation predicts that the helicity-independent recoil polarization for elastic 
electron-proton scattering vanishes. Assuming that the two-photon contribution is negligible, we used this reac- 
tion to measure the false instrumental asymmetries arising from misalignment, detector or tracking inefficiencies, 
variations of pathlengths in the analyzer, and other mechanisms. We express the detection probability in the form 

fh(8 'fpp, 0f pp ) = /o(6*fpp, <f>f PP )7r- (1 - he x sin cf) lpp + he y cos 4> lpp 

Ztt 

+ a Q sin (fopp + b cos cj> ipp + c sin 2cj> ipp + d cos 2cj> ipp ) (20) 

where the coefficients (ao, bo, Co, do) parametrize the false asymmetry while the coefficients (e Xl s y ) depend upon the 
helicity-dependent recoil polarization and the FPP analyzing power. Thus, the false asymmetry coefficients are 
obtained by Fourier analysis of 

^ + "t ^ = TT-C 1 + = 7T C 1 + a osin0 fp p + b cos(j) {pp + c o sin20 fpp + d cos 20 fpp ) (21) 
Z Jo Ztt Ztt 

Data for elastic scattering were taken at five proton momenta between 0.785 and 0.851 GeV/c and at 1.066 and 
1.188 GeV/c. The dependence of false asymmetries upon S are shown in Fig. ^5] for 5° < 0f pp < 20°. Only the 
coefficient of cos</>f pp shows a significant dependence upon S that can be attributed, in part, to its correlation with 
the average pathlcngth of scattered particles in the carbon analyzer. This dependence was parametrized by a linear 
function. The other three coefficients are essentially independent of 5 with average values less than 1%. No significant 
dependence upon proton central momentum is apparent over this range. 



F. Background subtraction 

The polarization response functions were corrected for two types of background: the elastic radiative tail and 
accidental coincidences. Corrections for the unresolved contribution of the elastic radiative tail were made using the 
likelihood function because those contributions varied strongly with both 8 and cf>. Thus, we generalize Eq. (|13fl to 

e a = A y (6 !pp ) ^T(/iSiX } + hS^Tf) ( 22 ) 



where f\ + /2 = 1 and where is the polarization for A excitation and T^ 1 is the polarization for elastic scattering. 
Note that the spin transformation matrix for elastic scattering differs from S^- 1 ' for pion production because 
the polarization vectors for the two reactions are normally evaluated in different frames. Thus, omits Rw (see 
Appendix and assumes that the proton emerges parallel to q. The relative weights depend upon the (W, Q 2 , x, (j>) 
bin and were obtained by fitting the missing-mass distributions for each bin with appreciable elastic contamination. 
This contribution is actually very small and is only visible for the 8 cm = —90° setting. Generalization of the likelihood 
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FIG. 13: (Color online) The dependence of measured Fourier coefficients for false asymmetry upon spectrometer relative 
momentum are shown with average values (do, Co, do) or a linear fit (bo) as blue lines. 



formula, Eq. I|18f) . is straightforward. The elastic polarizations were computed from the parametrizations of Ge p and 
Gmp found in Refs. |5fjl l57|. but the results are insensitive to the small differences between models. 

Accidental background was subtracted by analyzing both in-time and out-time events in the same manner. For 
polarization we obtain 

P = Pp \ rPb (23a) 

{spf _ w ; wp (23b) 

where P p ± SP p is the measurement for the in-time region, i 5 Pi, is the result for the out-time region, and r is the 
ratio between the widths of these regions. Similarly, for response functions we obtain 



R = Rp-rR b (24a) 
(SR) 2 = {6R P ) 2 + (rSR b ) 2 (24b) 



where R p ± SR P and Rb ± SRi, are obtained for in-time and out-time regions, respectively. The effect of background 
subtraction is generally difficult to discern in standard figures and is always much less than the statistical uncertainty 
in these measurements. 



G. Pseudodata tests 



The analysis procedures were tested using pseudodata. For each accepted event, response functions and polarizations 
at the target were computed based upon the MAID2000 model. The observed polar scattering angle 9f pp was retained 
but the azimuthal scattering angle tpp was sampled randomly. This value of </>f pp was retained if the likelihood L 
calculated according to Eq. Ijl2|l wa s greater than the next random number thrown and rejected otherwise. This 
procedure was iterated until a value of </) tpp was selected. Contributions to V and A were then accumulated and the 
pseudodata were analyzed in the same manner as real data. 

These tests demonstrate that model input for response functions is recovered within statistical uncertainties, but 
that there are sometimes inconsistencies in the <fi dependence of polarization data. This problem arises because 
relatively large bins in </> are needed to obtain useful statistical precision, but some of the spin-transport matrix 
elements can exhibit broad distributions with respect to other variables in part due to kinematic focusing in the 
lab frame. Under those conditions the acceptance-averaged polarization can differ appreciably from model values 
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for the central kinematics of a bin. These difficulties arc much smaller for response functions because binning with 
respect to cj) is not necessary; all 4> values contribute to the determination of a response functions and their coefficients 
are evaluated properly for each event. Explicit eventwise weighting with the leading factors of sin 6* also reduces 
the effects of acceptance averaging on the response functions as defined in Eq. ©. Therefore, we focus upon the 
response-function data and do not consider polarization binned with respect to <f> further. A more detailed report on 
the pseudodata analysis is provided in Ref. |5Sj . 



H. Acceptance averaging 

Multipole amplitudes and Legendre coefficients are functions of (W, Q 2 ), but the acceptance averaged (W, Q 2 ) 
depend upon x. Consequently extraction of these quantities from angular distributions can be distorted by the x 
variations of (W,Q 2 ). Such distortions can artificially enhance terms for large £„. Two methods for compensating 
for such distortions have been tested using both pseudodata and real data. The additive method is based upon the 
first-order expansion 

R(W, Q 2 ,x, e) = R(W, Q 2 ,x, e) - ^(W - W) - j^(Q* - Q 2 ) (25) 

where overlines indicate acceptance averaging and where the derivatives are evaluated at central kinematics using 
a model, such as MAID. For this experiment Q 2 — Q 2 tends to be much more important than W — W . Additive 
kinematical corrections have the advantage that variations of both W and Q 2 can be accommodated, but this procedure 
has the disadvantage that it relies upon a model and we have no model that provides a uniformly good fit to all of the 
response functions. While that is not a problem for pseudodata, the use of an inaccurate model to make kinematical 
corrections to real data could introduce more serious errors than it corrects. Therefore, a second procedure based 
upon form factors was tested. Assuming that all response functions share a common form factor, and that kinematical 
corrections are dominated by the x-dependence of Q 2 , we postulate 

R{W,Q 2 ,x) = R(W,Q^,x) [g(Q 2 )/G(Q^)Y (26) 

and approximate G(Q 2 ) by the usual dipole form factor G D (Q 2 J_= (1 + Q 2 /A 2 )~ 2 where A 2 = 0.71 (GeV/c) 2 . This 
multiplicative procedure does not compensate for variations of W, but for this experiment these variations are much 
smaller than those for Q 2 . 

Figure ITU compares pseudodata with multiplicative kinematical corrections with the model for central kinematics. 
The open squares show raw response functions extracted from pseudodata while open red circles show acceptance- 
averaged response functions from MAID2000. The agreement between these data sets, modulo statistical fluctuations, 
demonstrates the internal consistency of the simulation/analysis program. However, the x dependence of Q 2 can pro- 
duce significant systematic deviations from the input model (solid curves) evaluated at central kinematics, especially 
for R^ T - Recognizing that Q 2 0.94 for x > 0.5 or 1.06 (GeV/c) 2 for x < —0.5, we observe that the pseudodata for 
i?^ T do cluster around the model curve for the appropriate Q 2 . Even the abrupt transition across x = is reproduced. 
The solid circles show that "centered" pseudodata adjusted according to Eq. i|26|) cluster better around the model 
curves for central kinematics. Therefore, distortion of multipole amplitudes by the x dependence of Q 2 should be 
minimized by fitting centered data. We find that the multiplicative corrections move the pseudodata in the directions 
indicated by ratios between acceptance averaged calculations and those for central kinematics. The net effect is to 
reduce the scatter in the pseudodata and to remove many systematic deviations attributable to the x dependence 
of Q 2 . On the other hand, we have the qualitative impression that the corrections are sometimes a little too large, 
although no attempt has been made to quantify that impression. One could reduce the size of the multiplicative 
correction by increasing the dipole mass A, but the N — > A form factor is actually steeper (smaller A) than the 
standard dipole form factor. Furthermore, changes to the correction by replacing the dipole with a parameterization 
of the N — > A form factor are quite small. 

Therefore, we adopted the multiplicative correction based upon the dipole form factor as the standard method for 
bin centering. The figures in the remainder of this paper show recoil-polarization response functions plotted at x with 
bin-centering corrections for Q 2 . Legendre and multipole fits were made to the data in this form. 
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FIG. 14: (Color online) Pseudodata data for response functions at W = 1.23 ± 0.01 GeV and Q 2 = 1.0 ± 0.2 (GeV/c) 2 
are compared with the input model (MAID2000) at the central kinematics (solid curves) and with neighboring values of Q 2 
representative of the x dependence of acceptance averaging; dashed curves show Q 2 = 0.94 and dashed-dotted curves show 
Q 2 — 1.06 (GeV/c) 2 . Acceptance-averaged calculations are shown as red open circles and pseudodata as open squares. Filled 
circles show pseudodata with multiplicative kinematical corrections based upon the dipole form factor. 



I. Systematic uncertainties 



1. Response functions 



There are several types of normalization uncertainty that affect the response-function data. These include uncer- 
tainties in the unpolarized differential cross section used to normalize the likelihood, the FPP analyzing power, and for 
helicity-dependent responses the beam polarization. Although these systematic uncertainties do vary to some degree 
with spectrometer settings, beam conditions, and time, those variations are small compared with the statistical uncer- 
tainties. Therefore, we believe it sufficient to estimate average systematic uncertainties for those quantities without 
tracking the propagation of particular settings through event sorting. The typical systematic uncertainties in the 
differential cross section data are about ±3% point-to-point, so we assume that the uncertainty in the parametrized 
cross section used in the likelihood analysis is also about ±3%. Similarly, the systematic relative uncertainty for 
Compton measurements of beam polarization is estimated to be about 1.4% |4fJ. Finally, the relative uncertainties in 
average analyzing power reported by Punjabi el al. are in the 1 — 2% range. Because we do not consider thickness 
or momentum variations, we adopt a fairly conservative estimate of SA y /A y = 0.02. Therefore, the normalization 
uncertainties are approximately ±3.6% for helicity-independent or ±3.9% for helicity-dependent response functions. 

The evaluation of other types of systematic error requires replaying the data subject to a perturbation of one 
of the analysis parameters. Thus, the uncertainty due to subtraction of the elastic background was obtained by 
comparing replays with and without that subtraction. Because the contamination fractions binned in <j) were difficult 
to determine, we assumed their relative uncertainties to be 100% and estimated the corresponding uncertainties in 
response functions as 

SR a = \R^-R^\ (27) 

where Ra and Ra^ represent response function a with and without elastic subtraction. Similarly, the uncertainty in 
corrections for false asymmetry were estimated as 

SR a = 0.1\R^ -R^\ (28) 

where the relative uncertainty in false asymmetry was estimated to be ±10% and is multiplied by the difference in 
response functions obtained with and without false asymmetry in the likelihood function. 

A similar procedure was also applied for the spin transport matrix. The sensitivity of response functions to 
uncertainties in the spin rotation matrix is illustrated in Fig. ^|for W — 1.23 GeV. The open black points were 
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FIG. 15: (Color online) Sensitivity ol response functions for W = 1.23 GeV to uncertainties in the spin rotation matrix. 
Open black points were obtained using the COSY model while filled green points are based upon the Pentchev model. For 
the Pentchev analysis inner error bars are statistical while outer error bars include systematic errors due to uncertainties in 
precession angles and optical matrix elements; however, the systematic errors are generally too small to see. 

obtained using the COSY model while the filled green points were obtained using a simpler geometrical model 
by Pentchev in terms of 6 parameters consisting of two trajectory angles and four matrix elements coupling spin 
components [3?], EH- Obviously, this geometrical model accurately reproduces the COSY model. Therefore, we can 
estimate systematic errors in response functions due to uncertainties in the spin rotation by comparing results obtained 
from independent perturbations of each of the 6 parameters of the Pentchev model by its estimated uncertainty and 
combining in quadrature differences with respect to nominal parameters. We use the same systematic uncertainties 
for those parameters as in Refs. |37l I5f3j . The green error bars in Fig. 1151 include an inner statistical portion shown 
with endcaps and a total error without endcaps. However, rarely can one discern the systematic contribution to 
the total bar because the composite contribution of spin rotation errors is almost always small relative to statistical 
uncertainties. 

The net systematic uncertainties in response functions, consisting of the quadrature sum all contributions discussed 
in this section, tends to be dominated by the normalizations and is almost always small compared with statistical 
uncertainties. Many of the figures show statistical uncertainties as inner error bars with endcaps and total uncertainitcs 
as outer error bars without endcaps. The systematic contributions are occassionally visible for Legendre coefficients 
or multipole amplitudes, but are rarely visible for response functions. 



2. Legendre coefficients and multipole amplitudes 

Let y represent a fitted quantity, such as a quadrupole ratio, Legendre coefficient, or multipole amplitude and let 
(3 represent a calibration parameter or a scale factor applied to one of the corrections, such as false asymmetry. We 
estimate the systematic uncertainty (Sy) sys by adding in quadrature numerical derivatives 

(6y) 2 sys = £ (||%) 2 = £ Wi + Sfh) - y(A)) 2 (29) 

i. ' i 

estimated by performing a series of fits in which each calibration parameter is perturbed in turn. Therefore, Legendre 
and multipole fits to such data sets begin with the results of the best fit for nominal calibration parameters and 
usually require only a few iterations to determine small displacements of the minimum on the \ 2 hypersurface. We 
assume that the desired local minimum is related to the best fit by a small distortion of the % 2 hypersurface produced 
by small changes in the data set due to perturbation of an analysis parameter. By starting with the nominal best 
fit, we minimize the chance that the fitting procedure might find a different local minimum. With enough care in 
the fitting procedure, we find that changes in fitted Legendre coefficients or multipole amplitudes due to variation of 
spin-rotation parameters or omission of false asymmetry or elastic subtraction are typically small. 
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In addition to systematic uncertainties considered in the previous section, Legendre and multipole analyses also 
include an estimate of the uncertainty due to the kinematic or bin-centering correction. The customary dipole form 
factor should describe the Q 2 dependence of nonresonant contributions fairly well, but the N — > A form factor is 
known to have a more rapid Q 2 dependence. The best description is probably intermediate between these models. We 
estimated the uncertainty in fitted Legendre coeffcients and multipole amplitudes due to the choice of bin-centering 
form factor by comparing fits for the dipole and N — ► A form factors, assigning a systematic uncertainty equal to 
the difference between the two fits. The dipole form factor is given by Gd{Q 2 ) = (1 + Q 2 /A 2 )~ 2 where A 2 = 0.71 
(GeV/c) 2 . For the N — > A form factor we use the Sato-Lee parametrization Gna = (1 + aQ 2 ) exp (—bQ 2 )G d with 
a = 0.154 and b = 0.166 (GeV/c) -2 |6(j- However, the difference between these form factors over the range of Q 2 — Q 2 
for this experiment is too small to produce a visible difference in the projected data or fitted angular distributions. 

The systematic uncertainties in fitted Legendre coefficients and multipole amplitudes contain a total of 12 contribu- 
tions added in quadrature, each requiring a fit to the relevant data set. Variations of the cross section, FPP analyzing 
power, beam polarization, bin centering, false asymmetry, and elastic subtraction are all compared with the best fit 
for data obtained using COSY spin rotation. The six contributions to the spin rotation uncertainty are estimated 
using differences with respect to data based upon the nominal Pentchev model. The net systematic errors in these 
quantities are generally small compared with statistical uncertainties. Figures showing fitted quantities with statisti- 
cal and total errors bars can be found in the separate reports on Legendre coeffcients an d multipole amplitudes 

J. Summary of experimental data 

Near the middle of our (W, Q 2 ) acceptance, we have obtained complete angular distributions for 16 response 
functions, 14 separated plus 2 Rosenbluth combinations for e ~ 0.95. The angular coverage and statistical precision 
are clearly best in the central W range, 1.21 < W < 1.29 GcV. Data tables are on deposit with EPAPS |95J. These 
tables give both raw data and bin-centered data with both nominal and acceptance-averaged kinematics. Tables of 
Legendre coefficients and multipole amplitudes are included also. 

VI. RESULTS 
A. Comparison with Models 

In this section we compare our data for response functions with calculations using four recent models. We provide 
very brief summaries of the models and refer to original sources for more detailed information. A recent review of 
these and related models has also been provided by Burkert and Lee j6^| ■ 

The SAID model [l3L l64| parametrizes a photoexcitation multipole amplitude A in the form 

A = (A B + A Q ){1 + iUn) + A r Un + (C + iD)(lmt nN - % N ) (30) 

where t^N is a i-matrix fit to irN elastic scattering data that enforces the Fermi- Watson theorem [g^ below the two- 
pion threshold, Ar is parametrized as a polynomial in E n with the correct threshold behavior for each partial wave, 
Ab is a partial wave of the pseudoscalar Born amplitude, and Aq is parametrized using Legendre functions of the 
second kind. Recent extensions also include energy-dependent polynomials C and D. Electroexcitation amplitudes 
also include form factors in Q 2 . We are now using the WI03 version of SAID j^| . Multipole amplitudes were projected 
from helicity amplitudes using the formalism in Appendix [HI 

The Mainz unitary isobar model |6jfl , known as MAID, parametrizes resonant contributions to multipole amplitudes 
using the Breit-Wigner form 

F W 

A = A(Q 2 )C„ nM W) w2r _ ^ 2 * iWRT J *»(W) (31) 

where Wr is the resonance mass, r to t is its total width at resonance, C^n is an isospin factor, and A is a form factor. 
The W dependence of the electroexcitation vertex and its pseudothreshold behavior is represented by / 7 jv while f n N 
describes the R — > irN decay in terms of an energy-dependent partial width, r n pf(W), and appropriate phase-space 
and penetrability factors. Nonresonant amplitudes are computed using Born and vector- meson diagrams with a mixed 
itNN coupling that interpolates between pseudovector coupling at low cm momentum, p w , and pseudoscalar coupling 
at high Pit- Background amplitudes are unitarized with the (1 + U^n) factor, as above, while resonant contributions 
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are unitarized by adjusting the phase ip such that the total phase of the resonant contribution is given by the SAID 
partial-wave analysis for irN elastic scattering. Thus, ip depends upon both W and Q 2 and varies with multipole. 
The event generator used for data analysis employed MAID2000, but here we will also show calculations using the 
updated MAID2003 version [H HI . 

The Dubna-Mainz-Taipei (DMT) model [H E3 

is based upon MAID but employs a more sophisticated analysis 
of ttN rescattering. Whereas MAID employs a /T-matrix approximation for the background contribution to the 
i-matrix, DMT includes off-shell contributions in the form of a principal-value integral. Both models use similar 
Breit-Wigner parametrizations for resonances, but the electroexcitation amplitudes for MAID should be interpreted 
as "dressed" while for DMT the resonant amplitudes are considered "bare" because the ttN rescattering terms account 
for background contributions to resonant multipoles. 

The Sato-Lee (SL) model [60| is formulated in terms of energy-independent effective Hamiltonian and current 
operators. This dynamical model provides coupled equations for the ttN and 'jN reactions that automatically satisfy 
unitarity. The potential governing pion rescattering is optimized to reproduce ttN elastic scattering data. By means 
of a unitary transformation one can distinguish between the electroexcitation amplitudes for the N — > A transition 
and the contributions of the pion cloud and rescattering mechanisms. Although differing in detail, both the DMT 
and SL analyses conclude that the pion cloud is responsible for enhancing the Mi + amplitude relative to the quark 
model and for most of the observed quadrupole strength. Thus, these models suggest that the intrinsic deformation is 
rather small. Note that the SL model omits higher resonances and is limited to W ^ 1.4 GeV while the DMT model 
reaches larger W by including contributions of higher resonances based upon MAID2000. 

The data for response functions are compared in Figs. 1161251 with calculations based upon these models. The 
response functions in the first two columns, described as R-type, depend upon real parts of interference products 
while those in the last two columns, described as I-type, depend upon imaginary parts. Although the first three 
response functions in column 1 and the last in column 3 have been observed before, the other 12 response functions 
have been observed for the first time in this experiment. As a general rule we find that variations among the models are 
usually greater for I-type than for R-type response functions, although R TT also shows significant model dependence. 
When W « Ma, R-type responses arc largely determined by the relatively well-known multipole amplitudes for the A 
resonance while I-type responses require interference with nonresonant background or tails of nondominant resonances 
that are constrained less well by previous data. For both types the variations among models are typically smallest for 
W near and below Ma and increase with W above the A resonance. By the time we reach W ~ 1.3 GeV, variations 
among models become large even for R-type responses. Above the A resonance the magnitudes for many of the SL 
response functions decrease faster than the data as W increases, presumably due to neglect of higher resonances. 
Conversely, some of the DMT response become too strong as W increases, notably Rl+t, Rlt, and Rtt- SAID 
appears to be the least accurate of these models, especially for LT response functions. We speculate that this problem 
might be related to the use of pseudoscalar coupling for Born amplitudes. Among these models, MAID2003 seems to 
provide the best overall description of the data, but none provides a uniformly good fit. 

B. Legendre Analysis 

In the present representation, the response functions should be polynomials in x of relatively low order, especially 
if the assumption of M%+ dominance is valid near the A resonance. According to that assumption, one expects Rtt 
to be constant, Rlt, R'lt> R-l+t> Rtt' an< ^ Rtt t° be linear, and only R'l T , R'tt an( ^ Rlt t° be cubic; the others 
are expected to be quadratic. Indeed, at W — 1.23 GeV we find that the data for Rtt are almost constant and 
those for Rlt are almost linear, with deviations from these simple behaviors that are qualitatively consistent with the 
departures of the models from Mi + dominance. Similarly, at W = 1.23 GeV R' LT and R TT appear to be consistent 
with linear behavior despite the somewhat larger experimental uncertainties. However, model calculations for these 
responses show larger deviations with respect to M\+ dominance because imaginary parts of interference products are 
more sensitive to nonresonant mechanisms and tails of nondominant resonances. Finally, although R'[ T displays cubic 
behavior, R'^ T appears to be almost linear because the |Mi + | 2 contribution to its linear coefficient dominates the 
polynomial. However, these simple rules deteriorate rapidly as W increases and Rlt, Rl+t an d Rtt data develop 
strong curvatures. Furthermore, model calculations also show significant departures from these simple behaviors; for 
example, the curvature of Rtt calculations is often appreciable. 

Representative Legendre fits are compared in Figs. I26I3UI with data for response functions. A more complete 
set of figures and comparisons between fitted and predicted Legendre coefficients can be found in Ref. [fj]]]. Note 
that these fits employed the [x, <f) distribution for the differential cross section and beam analyzing power together 
with the data for recoil-polarization response functions. Thus, although Rl+t and Rtt could not be separated for 
x w directly, the Legendre fits to those response functions are determined well in this region nonetheless. The 
dashed curves are limited to the sp truncation while solid curves include additional terms in response functions 
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FIG. 16: (Color online) Data for response functions at W = 1.17 GeV are compared with selected models: MAID2003 (solid red), 
DMT (dashed green), SAID (dash-dotted blue), and SL (dotted cyan). Inner error bars with endcaps are statistical while outer 
error bars without endcaps include systematic uncertainties; however, the systematic contributions are often indistinguishable. 
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FIG. 17: (Color online) Data for response functions at W = 1.19 GeV are compared with selected models. See Fig. 1161 for 
legend. 



for which the sp fits appear to be systematically deficient over a range of W. Most of the response functions can 
be fit well with the truncation based upon M\+ dominance, but Rlt, Rtt, Rl+Tj an d Rtt generally require an 
extra term that reveals additional contributions. However, it is not immediately obvious whether those additional 
contributions arise from < 1 terms that do not involve M\+ or whether they require participation of higher 
partial waves. It is also important to recognize that even when Legendre expansions limited by Mi+ dominance 
do fit the data well, considerable violation of this assumption may still be present. Legendre fits are made to the 
data for each response function independently and ignore the correlations between response functions required by 
the expansions listed in Tables 1111111 More importantly, correct prediction of the number of significant Legendre 
coefficients using this truncation scheme does not ensure that the relationships between the values of those coefficients 
and the underlying multipolc amplitudes would be correct. A detailed study of the truncation errors in the Legendre 
analysis of unpolarized response functions and their consequences for simplistic extraction of multipolc amplitudes 
is provided in Ref. [7^ . Therefore, a more rigorous analysis that fits the multipole amplitudes direcly, without the 
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FIG. 18: 

legend. 



(Color online) Data for response functions at W = 1.21 GeV are compared with selected models. See Fig. 1161 for 
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FIG. 19: 
legend. 



(Color online) Data for response functions at W = 1.23 GeV are compared with selected models. See Fig. 1161 for 



mediation of Legendre coefficients, is presented in Sec. IVI CI 

(3/2) (3/2) 

Most previous extractions of R EM and R SM employed truncated Legendre analysis of the unpolarized differential 
cross section. To the extent that Mi+ dominance and sp truncation apply for W ~ Ma, one can define 
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FIG. 20: 
legend. 



(Color online) Data for response functions at W = 1.25 GeV are compared with selected models. See Fig. 1161 for 




FIG. 21: 

legend. 



(Color online) Data for response functions at W = 1.27 GeV are compared with selected models. See Fig. 1161 for 



are W-dependent combinations of Legendre coefficients for a particular charge state. To obtain the desired quantities 
for the A(1232) resonance, these quantities are evaluated at W = Ma and one needs to correct for the isospin-1/2 
contamination of the pir° channel. The isospin correction is expected to be small and has not been made. The 
reliability of Eq. JHSJl will be evaluated in Sec. IVf 1 Al 

Figure |^] compares values for ^ and R^m ^ obtained from Legendre fits to the differential cross section data 
for Q 2 = 1.0 (GeV/c) 2 with model calculations obtained from Eq. (|33|l where the Legendre coefficients were obtained 

by numerical integration. Although these quantities approximate R^m^ an< ^ ^SM oru y f° r W ~ their W 

dependence offers some insight into the model dependence of the traditional truncated Legendre expansion. The 
open circles were obtained using the M\ + truncation while the filled circles vary an additional term for each response 
function in order to improve the fits for larger W. Recall that Fig. ^2 demonstrates that the data are sensitive to at 
least one additional term per response function and that model calculations predict significant Legendre coefficients 
beyond M\ + dominance. Although the uncertainties increase for W > 1.3 GeV because the angular acceptance 
becomes too limited, we find that both analyses are qualitatively consistent for an appreciable range of W around 
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legend. 



(Color online) Data for response functions at W = 1.29 GeV are compared with selected models. See Fig. 1161 for 




FIG. 23: 
legend. 



(Color online) Data for response functions at W = 1.31 GeV are compared with selected models. See Fig. 1161 for 



A/a- The data for these quantities are relatively smooth with W dependencies that are similar to model calculations 
of the same quantities, whether or not these quantities are adequate approximations to the desired quadrupole ratios. 

The spread among models is smallest near the physical mass but remains appreciable for R^ M ^ , for which SAID 
differs significantly from both the data and the other models shown. Although the SL slope is somewhat too small 
compared with data near Ma, the other models provide a qualitatively consistent description of the W dependence 
of Rs m ■ In the central W region the experimental results are practically independent of truncation scheme and are 
in reasonable agreement with the MAID or DMT models, but the positive Rem values for SAID when W » Ma 
disagree sharply with the data. For larger W the SL calculation for Rsm is much flatter than the data, probably 
because higher resonances are omitted. Although it is more difficult to obtain unambiguous Rem results for W > 1.31 
GeV, data based upon the sp truncation remain in reasonable agreement with model calculations based upon the 
same truncation scheme. 

The results for W = 1.23 GeV are listed in TablelVland are practically independent of truncation scheme — the slight 
variation in Rem is within the estimated statistical uncertainty. We also list in Tablefvlthe values obtained by Joo et 
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FIG. 24: 
legend. 



(Color online) Data for response functions at W = 1.33 GeV are compared with selected models. See Fig. 1161 for 
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FIG. 25: 
legend. 



(Color online) Data for response functions at W = 1.35 GeV are compared with selected models. See Fig. 1161 for 



al. [lsf using the sp truncation, averaging with respect to neighboring Q 2 bins. Their results are consistent with ours 



for but are significantly larger for Rg M . Note that Joo et al. estimated that truncation errors in determination 

of quadrupole ratios were no more than 0.5% for SMR or 0.7% for EMR in absolute terms. While we agree that 
truncation of Legendre fits to the number of terms in the sp model has little effect upon fitted values for Rem or Rsm, 
the discussion in Sec. I Vlll demonstrates that the underlying assumptions of the traditional Legendre analysis do not 
provide adequate approximations to the quadrupole ratios at the present level of experimental precision. Therefore, 
the next section presents a more rigorous analysis based upon multipole fits that does not assume sp truncation or 
M\ + dominance. 

This analysis for the quadrupole ratios is based upon unpolarized cross sections and does not exploit any of the new 
recoil-polarization data. Examination of Tables IIIIIII shows that many other combinations of Legendre amplitudes 
could also provide ReEi + Mi + and ReSi + M^ + — if the sp truncation is valid these quantities should be highly 
overdetermined. Thus, one could obtain KeEi + Ali + and Re5i+M* + using a least-squares analysis of the entire set of 
Legendre coefficients for W = Ma and measure the reliability of the truncation scheme by x 2 - However, it is already 
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FIG. 26: Data for response functions at W = 1.21 GeV are compared with Legendre fits in the sp truncation (dashed) and 
with a few extra terms as needed (solid) . Inner error bars with endcaps are statistical; outer error bars without endcaps include 
systematic uncertainties. 






LT(1)' '- 






. ... i .... i . . 


TT(1)' ' ■ 




1 ' ' ' , . , . , | , , . . 


i i 

■ TT(t) 


i 


i*^ 





1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



FIG. 27: Data for response functions at W — 1.23 GeV are compared with Legendre fits. See Fig. 1261 for legend. 

clear from the Legendre fits that terms beyond sp truncation are needed for some of the response functions even for 
W » Ma- Furthermore, it is desirable to obtain the W dependencies of both the real and the imaginary parts of the 
multipole amplitudes. Therefore, we forgo further study of the consistency of sp truncation and proceed directly to 
a multipole analysis that exploits the new response functions. 

C. Multipole Analysis 



Let 



MW, Q 2 ) = A (0) (W, Q 2 ) + AA^W, Q 2 ) 



(34) 



represent either the real or the imaginary part of one of the multipole amplitudes (M^, Ey, or Sij) where is 
a baseline amplitude obtained from a suitable model while AAi is a variable to be fit to the data. To minimize 
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FIG. 28: Data for response functions at W — 1.25 GeV are compared with Legendre fits. See Fig. 1261 for legend. 
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FIG. 29: Data for response functions at W = 1.27 GeV are compared with Legendre fits. See Fig. |2§lfor legend. 



theoretical bias, our standard fits employ a baseline model consisting of Born terms for pseudovector ttNN coupling 
plus p and ui exchange; see Appendix lEl for details. To test the sensitivity of fitted multipolc amplitudes to neglect 
of tails of nondominant resonances or to variations of nonresonant contributions, we have also performed fits using 
MAID2003, DMT, SL, or SAID as baseline models. Note that some of the AAi parameters are relatively large 
when using the Born baseline that contains no information about the A(1232) resonance while the fitting parameters 
for more sophisticated baseline models represent small corrections to the specified model. Nevertheless, we have 
demonstrated that fitted multipole amplitudes are practically independent of the baseline model; see Ref. for 
details and figures. Both Legendre and multipole analyses were performed using the EPIPRDD program [73| . 

Fits were performed for each W bin to data consisting of the (x, <j>) distributions of differential cross section and 
beam analyzing power plus x distributions for 12 recoil-polarization response functions simultaneously. Fits using 
Born amplitudes for pseudovector coupling as the baseline model are shown in Figs. 13 214 II Several fits were performed 
to determine the maximum number of parameters that can be extracted without flattening the x 2 hypersurfaces too 
severely or encountering uncontrollable correlations among parameters. Dashed curves, designated sp, show fits that 
adjust real and imaginary parts of all s- and p-wave multipole amplitudes with higher partial waves constrained by 
the baseline model, here just real Born amplitudes without resonances. The fits designated spd3> also vary the real 
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FIG. 30: Data for response functions at W — 1.29 GeV are compared with Legendre fits. See Fig. 1261 for legend. 



TABLE V: Quadrupole ratios for W = 1.23 GeV at Q 2 — 1.0 (GeV/c) 2 . The reduced chisquare for the entire data set is labeled 
xl while the chisquare per point for differential cross section data is labeled Xjv(cr). 



method 


SMR, % 


EMR, % 


xl 


xU°) 


sp 


-6.07 ±0.11 


-2.04 ±0.13 


1.7 


1.6 


sp+ 


-6.11 ± 0.11 


-1.92 ±0.14 


1.5 


1.4 


Joo et al. a 


-7.4 ±0.4 


-1.8 ±0.4 







"Weighted average of Q 2 = 0.9 (GeV/c) 2 results for = 1.645 and 2.445 GeV from 0. 



parts of 2— multipoles and are shown as blue dotted curves. Fits designated spd vary real and imaginary parts of all 
amplitudes with £ v < 2 and are shown as green dash-dotted curves. Finally, the fits designated "final" are similar 
to the spd3 fits except that ImMi_ is held at baseline, here zero, for reasons discussed below. There are 14 free 
parameters for each W in an sp fit, 26 for an spd fit, 17 for an spd3 fit, or 16 for the final fit. For comparison, 
Legendre fits vary 50 free parameters in the central W region. 

These figures show that fitting just the s- and p-wave amplitudes, with a Born background, is already sufficient to 
obtain a good description of the data. Although fitting d-wave amplitudes, or a subset thereof, sometimes provides 
modest reductions of x 2 , there is little systematic evidence that modification of d-wave or higher multipoles is really 
necessary. However, it is also clear that variation of all £ n < 2 amplitudes offers too much freedom — the oscillations 
in green dash-dotted curves for W > 1.29 GeV or W < 1.19 GeV are not needed to fit the data, are implausible in 
amplitude, and change too much from one W to the next. 

The fitted multipole amplitudes are compared in Figs. I42l46l with several recent models 0, [n], [^J . In addition, 
the baseline Born amplitudes are shown by solid curves. Note that all multipole amplitudes are real in this model 
and there are no resonances; therefore, the starting conditions are quite poor and large adjustments to the initial 
parameters are required to fit the data. Nevertheless, the fits describe the data well, the fitted parameters generally 
display smooth W dependencies, and the characteristic resonance profiles emerge in the 1+ multipole amplitudes 
without coaching. Note that the sp and final fits began with baseline amplitudes but to improve stability the fits 
with more freedom began with the results of the final analysis. The uncertainties increase with the number of free 
parameters because the data do not adequately constrain multipoles for £ v > 2. Thus, we reject the full spd analysis 
because the uncertainties in its parameters are large and the resultant oscillations in calculated response functions are 
not warranted by the available data. The amplitudes and resulting fits for the other three analyses tend to be very 
similar except that there is a rather strong correlation between ImMi_ and Im<Si_ for small W that results in fitted 
values of opposite sign that are substantially larger than model predictions for W < 1.21 GeV. This correlation also 
appears to affect imaginary parts of 1+ amplitudes for W = 1.17 GeV. Evidently, the data for low W do not distinguish 
between ImMi_ and lmS±- well enough to fit both simultaneously. Therefore, our final analysis eliminates ImMi_ 
because models tend to predict stronger ImSi_ amplitudes and the sp fit also produces rather small ImAf^ values. 
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FIG. 31: (Color online) The truncated Legendre analysis for EMR and SMR at Q 2 = 1.0 (GeV/c) 2 is compared with MAID2003 
(solid red), DMT (dashed green), SAID (dash-dot blue), and SL (dotted cyan). Both theory and experiment employ Eq. I33H 
and approximate EMR and SMR at W = Ma, indicated by the vertical line. Open circles are fit according to sp truncation 
while filled circles permit additional freedom in the Legendre analysis. For filled circles, inner error bars with endcaps are 
statistical while outer error bars without endcaps include systematic uncertainties. 
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FIG. 32: (Color online) Multipole fits for W = 1.17 GeV at Q 2 = 1.0 (GeV/c) 2 using a Born baseline model. Dashed curves 
fit corrections to all s- and p-wave amplitudes, blue dotted curves also fit real parts of 2— multipoles, and green dash-dotted 
curves fit all s-, p-, and d-wave amplitudes. The solid red curves, considered the final fit, are similar to the blue dotted curves 
except that ImMi_ is absent. See text for further details. 



The uncertainties in fitted multipoles is reduced and the W dependencies are improved, especially for imaginary 1+ 
amplitudes, upon elimination of this redundant parameter. Furthermore, Figs. l32l41l demonstrate that elimination of 
ImAfi_ does not visibly reduce the quality of the fits to the data. The "final" analysis varies both real and imaginary 
parts of Si- and all 0+ and 1+ multipoles plus the real parts of M%- and all 2— multipoles for a total of 16 free 
parameters for each (W, Q 2 ) bin. 

There is rather little spread among models for Mi + amplitudes across the A(1232) resonance and our experimental 
amplitudes agree well with model calculations even when the fit is based upon a Born baseline model without reso- 
nances. The variation among models is also relatively small for 5i+ amplitudes and good agreement is obtained with 
data for ImS'i + , but for the real part the present data exhibit a steeper slope on the large W side. MAID2003, DMT, 
and SL calculations are similar for E\ + , but the current SAID calculations are substantially different and disagree 
with the data. Our results for Re£i+ agree relatively well with MAID2003, DMT, or SL but the Im£i + data do not 
show the node near W ~ 1.27 GeV predicted by those models; there is no sign change for W < 1.35 GeV. 

Among the models considered, MAID2003 tends to provide the best description of the recoil-polarization data, but 
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FIG. 33: 
legend. 



(Color online) Multipole fits for W = 1.19 GeV at Q 2 = 1.0 (GeV/c) 2 using a Born baseline model. See Fig 
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FIG. 34: 
legend. 



(Color online) Multipole fits for W = 1.21 GeV at Q 2 = 1.0 (GeV/c) 2 using a Born baseline model. See Fig. 1321 for 



it does not reproduce the R e LT , I^ LT , or i?£ T angular distributions (see Figs. I18I22|1 . These difficulties appear to 
arise primarily from the So + amplitudes. Whereas MAID2003 suggests a nearly constant ReSo+ amplitude in this W 
range, we find less negative results that are consistent with the negative slope in W suggested by SAID. Although 
the SL calculation crosses the ReSo+ data near the middle of the W range, it has the opposite slope. We also find a 
rather steep slope for IrjiiSo+- It is interesting to observe that the SAID model agrees best with the 0+ amplitudes 
even though it is worst, among these models, for Ei+, Rel— , and 52-. All of the models agree pretty well with the 
HcEq + data, but none reproduce the W dependence seen for Imi?o+- 

It is interesting to observe that there is a wide spread among the models for ReMi_ but that the data are closest 
to the Born model that omits the Roper resonance, which suggests that the transverse amplitude p Ai/ 2 is small. 
On the other hand, the fitted Re5i_ does differ from the Born model and suggests that there is a nonnegligible 
longitudinal contribution from the Roper consistent with a radial excitation. These models agree fairly well with the 
ImS*!- data but, with the exception of the large SAID prediction, are smaller than the Re5i_ data. Thus, it appears 
that excitation of the Roper resonance is primarily longitudinal at Q 2 — 1 (GeV/c) 2 . We will return to this issue in 
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FIG. 35: 
legend. 



(Color online) Multipole fits for W = 1.23 GeV at Q 2 = 1.0 (GeV/c) 2 using a Born baseline model. See Fig 
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FIG. 36: 
legend. 



(Color online) Multipole fits for W = 1.25 GeV at Q 2 = 1.0 (GeV/c) 2 using a Born baseline model. See Fig. 1321 for 



Sec. IVIIII using a unitary isobar model. 

The real 2— amplitudes are small in this W range, but the slope for Rci?2- appears to be determined well by these 
data and is in good agreement with Born, MAID2003, or DMT predictions. However, the predictions of the SAID 
model are much larger than the data for ReS^-. SAID also predicts significant oscillations in 2+ amplitudes that are 
absent in other models. Although we cannot fit the 2+ amplitudes accurately, the large = 2 amplitudes for SAID 
produce oscillations in many of the response functions that are not warranted by the data. 
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FIG. 37: 
legend. 



(Color online) Multipole fits for W = 1.27 GeV at Q 2 = 1.0 (GeV/c) 2 using a Born baseline model. See Fig 
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FIG. 38: 
legend. 



(Color online) Multipole fits for W = 1.29 GeV at Q 2 = 1.0 (GeV/c) 2 using a Born baseline model. See Fig. 1521 for 



D. Quadrupole ratios 

The quadrupole deformation parameters can now be obtained directly from the fitted multipole amplitudes using 

Ei+ Re£a+M 1 * + _ Re^i + ReMi + + ImEi + ImMi + 



Re 
Re 



Mi+ 
Si+ 
Mi+ 



1 

ReSi+M 



\M- 



i+ 



ReMi+ReMi + + ImMi+ImMi + 
RcS'i + RcAf 1+ + Im5i + ImMi+ 
ReMi+ReMi + + ImMi+ImMi + 



(35a) 
(35b) 



where the multipole analysis provides the real and imaginary parts of each amplitude separately. The W dependencies 
of quadrupole ratios for the pir° channel are shown in Fig. 1471 and the results at W = 1.23 GeV are listed in Table 
fVll the correction for the small isospin 1/2 contamination is discussed in Sec. IVIII The truncation dependence of the 

experimental results is relatively small for -Rg^ ' and most of the models are in good agreement with the data for 
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FIG. 39: 
legend. 



(Color online) Multipole fits for W = 1.31 GeV at Q 2 = 1.0 (GeV/c) 2 using a Born baseline model. See Fig 
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FIG. 40: 
legend. 



(Color online) Multipole fits for W = 1.33 GeV at Q 2 = 1.0 (GeV/c) 2 using a Born baseline model. See Fig. 1321 for 



W ~ Ma, but SMR is significantly stronger for SAID. Although the truncation dependence is larger for R%J^ data, 
the value at Ma still appears to be determined relatively well and is consistent with all of the models except SAID, 
which gives a much smaller value and at larger Q 2 sign opposite other models. The model calculations spread more 
rapidly for the electric than for the scalar ratio as the distance from Ma increases. 

Note, that if one defines Ma as the W where ReM-^ 2 -* = 0, then the quadrupole formulas in Eq. (|35|l reduce to 



D (3/2) T ,-,! 

K SM = lm5 1+ 



(3/2) 



(3/2) /ImM{f» 



(36a) 
(36b) 



for the isospin-3/2 channel. However, these formulas are unsuitable for data analysis because comparable nir + data 
are not available for isospin decomposition and because the appropriate value of Ma is not known precisely or 
independently of models. It would also be necessary, in principle, to interpolate the multipole data with respect to 
W. We employ Eq. (|35|l because it is independent of W, applies equally well to pir° or isospin-3/2, and does not 
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FIG. 41: (Color online) Multipole fits for W 
legend. 



1.35 GeV at Q 2 — 1.0 (GeV/c) 2 using a Born baseline model. See Fig. 
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FIG. 42: (Color online) Fitted 1+ multipole amplitudes using Born baseline and adjusting all l w amplitudes except ImMi_ plus 
real parts of 2— amplitudes. Inner error bars with endcaps are statistical; outer error bars without endcaps include systematic 
uncertainties. The baseline amplitudes are shown as black curves. Several other recent models are shown also: MAID2003 (red 
dashed), DMT (green dot-dashed), SAID (blue dotted), and SL (cyan short-dashed). 
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FIG. 43: (Color online) Fitted 0+ multipole amplitudes using Born baseline. See Fig. 1421 for legend. 
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FIG. 44: (Color online) Fitted 1— multipole amplitudes using Born baseline. See Fig. 1421 for legend. 
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FIG. 45: (Color online) Fitted 2— multipole amplitudes using Born baseline. See Fig. ^Jfor legend. 



require any model-dependent assumptions about Ma. Furthermore, because the energy dependence in Fig. 1471 is 
quite mild, Table IvTl simply lists values for the bin closest to Ma, namely W = 1.23 GeV. Small corrections for the 
energy dependence of these quantites are evaluated in Sec. IVII CI 

Table fWI evaluates the sensitivity of quadrupole ratios to the selection of adjustable multipoles. The uncertainties 
increase when higher partial waves that are not well-constrained by the data are permitted to vary. Above we argued 
that the best compromise is obtained by varying 0+, 1+, 1—, and real parts of 2— multipole amplitudes with 2+ and 
higher partial waves constrained by the baseline model. Elimination of ImMi_ further reduces the uncertainties in 
SMR, without affecting the quality of the fit, by suppressing its unresolvable correlation with ImSi_. As previously 
argued, we believe that elimination of ImAfi_ is justified by the prediction, by all models considered, that it is 
negligible in this energy range. The results in the last two lines of Table IvTl are practically identical for SMR, though 
with reduced uncertainty in the final line, while the change in EMR is within the estimated uncertainites. 
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FIG. 46: (Color online) Fitted 2+ multipole amplitudes using Born baseline. The open circles show results for the spd analysis. 
Only statistical errors are shown. See Fig. 1421 for curves. 
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FIG. 47: (Color online) Multipole analyses for EMR and SMR at Q 2 = 1.0 (GeV/c) 2 are compared with MAID2003 (solid 
red), DMT (dashed green), SAID (dash-dot blue), and SL (dotted cyan). The vertical line shows the physical mass, W = Ma- 
Open circles adjust £„ < 1 multipoles while filled circles represent the final fits. For filled circles, inner error bars with endcaps 
are statistical while outer error bars without endcaps include systematic uncertainties. 



TABLE VI: Quadrupole ratios for W = 1.23 GeV at Q 2 — 1.0 (GeV/c) 2 using the pseudovector Born baseline model. Only 
statistical uncertainties from fitting are given. 



variables 


SMR, % 


EMR, % 


xl 


o+,i+,i- 


-6.73 ±0.24 


-2.43 ±0.19 


1.69 


0±,l±,l-,2-,2± 


-6.95 ± 0.49 


-3.19 ±0.79 


1.64 


0+, 1 + , l-,Re2- 


-6.85 ± 0.27 


-2.73 ±0.20 


1.65 


above except ImMi_ 


-6.84 ±0.15 


-2.91 ± 0.19 


1.65 



E. Sensitivity to baseline model 



As mentioned above, the multipole fits are rather insensitive to the choice of baseline model. To illustrate this, 
Fig. 1481 compares fits to the response functions for W = 1.23 GeV based upon several baseline models; figures for 
other W bins are available in Ref. |62j |. The fits based upon Born, MAID2003, DMT, or SL models are practically 
indistinguishable. The fits based upon SAID display a more oscillatory structure that is not supported by the data in 
the middle of the W range where the precision is best. The oscillations are presumably due to relatively large Rei?2+ 
and Rc^2+ amplitudes that are not ameliorated by the current truncation scheme. Therefore, the data clearly require 
smaller 2± amplitudes than predicted by the SAID model. 

Similarly, the fitted multipole amplitudes are also rather insensitive to the choice of baseline model. Even the fits 
based upon SAID, starting from rather different initial conditions and with significantly larger fixed 2± amplitudes, 
converge upon essentially the same final results. For example, Table IVm lists the quadrupole ratios for W = 1.23 GeV 
based upon several choices of baseline models and using the "final" parameter space. All of the results are consistent 
except those using the SAID baseline, for which SMR is substantially higher and EMR lower than for other baseline 
models. However, the quality of the fit is also noticeably inferior even though the differences in \t are not impressive. 
Therefore, we conclude that this version of the SAID model does not provide a suitable baseline for multipole analysis 
and judge the sensitivity to uncertainties in the baseline model to be similar to the tabulated fitting uncertainties. 



TABLE VII: Dependence of quadrupole ratios for W = 1.23 GeV at Q 2 = 1.0 (GeV/c) 2 upon baseline model. All s- and p-wave 
amplitudes, except ImMi_, plus real 2— amplitudes were fit with other amplitudes given by the specified baseline model. 

baseline SMR, % EMR, % X? 

Born -6.84 ±0.15 -2.91 ±0.19 1.65 

MAID2003 -6.90 ±0.15 -2.79 ±0.19 1.67 

DMT -6.82 ±0.15 -2.70 ±0.19 1.67 

SL -6.79 ±0.15 -2.81 ±0.19 1.64 

SAID -7.38 ±0.15 -2.53 ±0.20 1.85 
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X X X X 

FIG. 48: (Color online) Sensitivity of multipole fits for W = 1.23 GeV at Q 2 = 1.0 (GeV/c) 2 to the choice of baseline model: 
Born (black solid), MAID2003 (red dashed), DMT (green dash-dotted), SAID (blue dotted), or SL (cyan short-dotted). 

VII. DISCUSSION 

A. Reliability of traditional Legendre analysis 

Both the Legendre and multipole analyses fit the data well, but they yield significantly different estimates of 
the N — ► A quadrupole ratios. The Legendre results are listed in the first line of Table I Villi and those for the 
multipole analysis in the second and fifth columns of the second line. Subsequent lines show model calculations for 
quadrupole ratios based upon several truncation schemes. The second and fifth columns are the proper ratios of 
multipole amplitudes while the remaining columns use the traditional estimators given by Eq. i|33|) with Legendre 
coefficients that were computed by numerical integration of response functions obtained from the indicated truncation 
of the multipole amplitudes with respect to £ n . Thus, the third and sixth columns represent the sp truncation while 
the fourth and seventh columns are practically complete with respect to £^. We placed the experimental Legendre 
results in the £ v < 5 columns because truncation is not possible experimentally. The model Legendre coefficients 
were computed without using Mi+ dominance but the corresponding traditional quadrupole estimators, Rsm and 
Rem, employ combinations that were derived under that assumption. The values of Rsm and Rem for £ n < 5 
obtained from the fitted multipole amplitudes are similar to those obtained from the fitted Legendre coefficients but 
are distinctly smaller than the fitted values for Rsm and Rem even though the fits to the cross section data are 
practically identical. The differences between £ n < 1 and £ n < 5 model calculations demonstrate that sp truncation 
is often a poor approximation to Rsm and Rem, especially for the latter. Although the sp truncation of Rsm is 
reasonable for the SAID and SL models and for the present multipole fit, it is inaccurate for the MAID2003 and DMT 
models. However, sp truncation of Rem is quite poor for all models considered. Furthermore, the correspondence 
between the traditional Legendre estimators and the actual quadrupole ratios also depends upon the requirement 
that Mi + appears in every term of the multipole expansion of Legendre coefficients. The differences between the 
£ n < 5 results and the actual quadrupole ratios demonstrate that the assumption of Mi+ dominance is not sufficiently 
accurate either. 

A more detailed study of truncation errors in the traditional Legendre analysis of N — > A quadrupole ratios has been 
provided in Ref. [T^- Truncation errors are especially severe for Rem where the contribution of ReMi_.E* + alone is 
approximately —40% of the leading term using MAID2003 pir° multipoles for our kinematics. Many other neglected 
terms are significant and the convergence is slow and model dependent. Furthermore, the contributions of Aq and A\ 
to Eq. (|33|l are not negligible, as assumed using M\ + dominance. The contribution of A\ can, in fact, have a large 
effect upon delicate cancellations within the numerator of Rem- Thus, Rosenbluth separation should be performed 
before using the Legendre method, especially for Rem, but none of the recent Legendre analyses [13, H, W& 12JJ have 
done so, including the present experiment. The convergence of Rsm is better, but its relative accuracy as an estimate 
of Rsm is still no better than about 20% [T^- Therefore, although the details are model dependent, it is clear that 
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TABLE VIII: Calculated quadrupole ratios for W = 1.23 GeV at Q 2 = 1.0 (GeV/c) 2 . Columns with ranges of £ v are based 
upon traditional estimators given by Eq. £2). 



model 


U SM > /o 


n SM , /o 
£r < 1 


rt SM , /o 

tr < 5 


EM > /o 


nip™ ) 0/ 
EM ' /0 

&r < 1 


p(pT°) (V 
"•EM ' /0 

< 5 


Legendre fit 






-6.11 






-1.92 


multipole fit 


-6.84 


-6.46 


-6.00 


-2.91 


-1.54 


-2.18 


MAID2003 


-6.73 


-6.37 


-5.63 


-1.65 


-0.57 


-1.12 


DMT 


-7.21 


-6.77 


-6.10 


-1.77 


-0.70 


-1.47 


SAID 


-8.71 


-7.78 


-7.88 


+0.17 


+1.96 


+0.22 


SL 


-6.59 


-6.69 


-6.58 


-2.29 


-1.29 


-1.58 



neither assumption of the traditional Legendre analysis is sufficiently accurate at the present levels of experimental 
precision and completeness. 

B. Isospin-1/2 contamination of EMR and SMR 

Separation of the isospin 1/2 and 3/2 contributions to the multipole amplitudes would require comparable data 
for the nn + channel, including angular distributions for either recoil or target polarization, which are presently 
unavailable. Fortunately, the isospin 1/2 contamination is expected to have relatively little effect upon the deter- 
mination of isospin 3/2 quadrupole ratios. According to MAID2003, one expects (Rsm,Rem) = (— 6.71%, —1.62%) 
at (W,Q 2 ) = (1.23,1.0) for isospin 3/2 compared with (-6.73%, -1.65%) for the pn° channel H|. Similarly, the 
quadrupole ratios for the pure N — > A contribution become (—6.73%, —1.53%) in the absence of background. Finally, 
if one attributes the correction terms for 1+ multipoles in Eq. (|34|l entirely to the A resonance, assuming that the 
Born baseline model is accurate, we would estimate 

6.81% (37a) 

3.12% (37b) 

in good agreement with the full results for the pir channel that include background. However, the fact that changes 
in Rem due to neglect of background are opposite for MAID2003 and the experimental multipole fit suggests that 
part of the fitted AEi + should probably be attributed to the background in the baseline model. Nevertheless, it 
appears that corrections for isospin 1/2 contamination are probably smaller than the present error bars. Therefore, 
this model-dependent correction has not been made. 

C. W dependence of EMR and SMR 

The W dependencies of quadrupole ratios obtained using both Legendre and multipole analyses are compared in 
Fig. 0^] with parabolic fits of the form 

2 

y = 5> fc G^-Af A ) fe (38) 

fc=0 

using a nominal value of Ma = 1.232 GeV. The fits were confined to the central region, indicated by dotted vertical 
lines, where this simple parametrization should suffice for interpolation. Both fits describe the data for the central 
region well. It is important to remember that the quadrupole ratios for the two types of analysis are different quantities 
and need not have the same shapes — the Legendre estimators are affected by e and by all partial waves while those 
from multipole analysis are not. The parabolic fits appear to extrapolate better for the multipole analysis than for 
the Legendre analysis but should still only be used in the central region. 

The expansion coefficients fitted using the weighted linear least-squares method are given in Table H3 where the 
data for quadrupole ratios are expressed in percent and where the multipole amplitudes are based upon the Born 
baseline. The eto parameters represent the best estimates of the quadrupole ratios at Ma but fits with only 2 degrees 
of freedom do not necessarily provide realistic uncertainties. Instead, we quote the largest of 5/2(1.23) , Sao, and 
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TABLE IX: Power series for quadmpole ratios in terms of (W — Ma) in units of GeV. The best estimate of the quadrupole 
ratios at W = Ma — 1.232 GeV is given for each analysis method by oo. 







EMR 






SMR 




method 


ao 


ai 


0,2 


ao 


ai 




Legendre 


-1.76 ± 0.19 


-3.77 


-793 


-5.87 ± 0.20 


53.2 


618 


multipole 


-2.87 ± 0.19 


-13.3 


-450 


-6.61 ± 0.18 


39.3 


749 



fi a oyxi where 5R(1.23) is the uncertainty in the single-energy fit to data for W = 1.23 GeV, Sao is the uncertainty in 
the value of ao according to the linear-least squares fit to the energy dependence, and % 2 is the reduced chisquare for 
that fit in the central region. The second-order terms are negligible between Ma and the nearest W bin and changes 
due to the linear terms are less than one standard deviation. For example, the fitted EMR and SMR values for the 
multipole analysis are —2.85 and —6.69% at W = 1.23 GeV. The remaining differences between these values and those 
listed on the last line of Table IVT1 are less than one standard deviation. Although the quadrupole ratios in Tabic 17x1 
are slightly smaller, they are consistent with those we reported in Ref. [3^|. The analysis in Ref. [3^ considered only 
a single energy, W — 1.23 GeV, while the present analysis fits the energy dependence within the central region. The 
small differences are partly due to the change in W from 1.23 to 1.232 GeV and partly due to statistical fluctuations 
of the data for 1.23 GeV relative to the average trends represented by the curves in Fig. 0^1 Therefore, we consider 
the interpolated values in Table llXl to be the best estimates of the quadrupole ratios for W — Ma- 



D. Relationship between Ge„ and Rsm 

Buchmann jT^ has derived a relationship 

_ qMN_ G En . v 

between Rsm for the N — > A transition and the neutron electric and magnetic Sachs form factors, Geu and Gun- 
Here q is virtual photon momentum in the cm frame. Deviations from this relationship were attributed to three-quark 
and higher-order currents and were estimated to be at the level of 1/N 2 , or about 10%. Figure [5U| compares this 
prediction with recent data 0, 0, 0, 0] where the band is based upon fitted neutron form factors from Ref. . 
Note that the growth of the band for Q 2 > 1.5 (GeV/c) 2 where Geu data are presently unavailable, is artificially 
limited by the use of a model with only two parameters. The Buchmann formula underestimates most of the Rsm 
data. The discrepancy of about 15% at Q 2 — 1 (GeV/c) 2 is similar to the estimated theoretical uncertainty, but this 
model predicts a nearly constant quadrupole ratio for larger Q 2 while the data show a steep slope. Note that for 
the high-Q 2 data we chose the effective Lagrangian analysis instead of the Legendre analysis from Ref. because 
truncation errors in the Legendre method are expected to increase with Q 2 j23 and the effective Lagrangian results 
are consistent with the MAID and DMT analyses of the same data [7(j . Although the Rsm slope is described well by 
dynamical models of 7rA^ rescattering, perturbative QCD predicts that Rsm should become constant asymptotically. 
Therefore, it would be of interest to extend measurements of Geu to higher Q 2 and to use model- independent multipole 
analysis of new polarization data for pion electroproduction to verify the apparent slope in Rsm- 
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FIG. 50: Comparison between Rsm data and Buchmann's formula (band) using fitted neutron form factors. Data: open square 
[T3 |. filled triangle |T^ |. open triangle [2l| . open circles [Ts| . crosses [Tfl| |: the filled circle is the present result. Small horizontal 
displacements are used to reduce clutter. Error bars include statistical and systematic but not model uncertainties. 



E. Sensitivity of Legendre coefficients to specific multipole amplitudes 

Even though sp truncation and M\ + dominance are not sufficiently accurate for quantitative analysis of the 
quadrupole ratios, that truncation can still provide qualitative insight into the sensitivity of selected Legendre coeffi- 
cients to particular multipole amplitudes. For example, Table ITTll shows us that Aq TI w SlmM-^Mi- while Fig. ^] 
shows that most models predict that Mi_ is nearly real and slowly varying over this range of W. Consequently, one 
expects the W dependence of A^ Tt to strongly resemble ImA/i+ and its amplitude to be proportional to ReMi_ and 
opposite in sign. Figure ISTl shows that these expectations are realized by the Legendre fit. The observation that the 
Ag Tt data are smaller than SAID, larger than DMT and SL, and in good agreement with MAID2003 predictions is 
consistent with the same pattern seen in Fig. 1441 for ReMi_ and with response function figures in Sec. IVI Al A similar 
correspondence is also observed for A^ +Tn , but the truncation is not as reliable because Rosenbluth separation is not 
available and model calculations have greater shape differences. The ReAfi_ amplitude also appears in A\ Tn but is 
again diluted. Therefore, the best sensitivity to ReMi_ is provided by the Rtp T response function. 

Similarly, model calculations suggest that S\- also varies relatively slowly and is nearly real in this W range, 
although neither feature is quite as accurate as for Mi_. Inspection of Table ITTll then suggests that the best sensitivity 
to Re«Si_ is offered by the -R£ T response function through its Ag Tn Legendre coefficient. The shape differences show 
that Mi + dominance is not as accurate for this Legendre coefficient, but Fig. shows that its W dependence does 
resemble that of ImMi + nonetheless. Thus, we find that the SAID prediction for Aq Tti is considerably too strong 
while those of MAID2003, DMT, and SL are too weak and that the same pattern is observed in Fig. ^]for ReSi_. 

Within the truncated Legendre expansion, Rei?o+ is isolated by Aq Tti , Ag Te , or Aj Tt , which should be equal 
modulo signs if Eq + were real and Mi+ dominance accurate. The A\ Tt coefficient is included in Fig. |^but the other 
figures are omitted and can be found in Ref. |61| . We do observe the expected pattern of signs and all are similar in 
shape to ImA/i + but their magnitudes do not conform to these simplistic predictions. Nevertheless, the relationships 
between data and model calculations for the Legendre coefficients are similar to those for Rei?o+ • 

The most complicated situation is Sq+ because both model predictions and fitted amplitudes show important imag- 
inary contributions to this nonresonant partial wave. Hence, measurement of «So+ amplitudes requires LT response 
functions of both R-type and I-type. The So+ contributions to the truncated Legendre expansion are isolated by 5 
R-type and 5 I-type coefficients but, although each group displays a relatively uniform shape with respect to W, there 
are significant differences in detail that show that the simple truncation is not especially accurate for these coefficients. 
Nevertheless, the fitted multipole amplitudes provide good fits to all of the LT response functions simultaneously. 
Therefore, polarization data provide the phase of Sq+ . 

Schmiede n 1221 speculated that the disagreement between SMR values for Q 2 ~ 0.13 (GeV/c) 2 obtaine d by Kalle- 
icher et al. |2l| using the forward pions versus those obtained using cross sections |l5f or polarizations |l4l |23| for 
forward protons might be explained using Re(So+/Mi + ) ps —0.14. However, recent models predict smaller positive 
values with relatively slow Q 2 dependence that are similar to —Rsm- Although it might appear that one could 
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FIG. 51: (Color online) Fitted Legendre coefficients for _R^ T are compared with MAID2003 (red solid), DMT (green dashed), 
SAID (blue dash-dotted), and SL (cyan dotted). In the simplest approximation, the first and second panels are — 3ReMi_ImMi+ 
and 3ReI?o+ImMi+, respectively. Inner error bars with endcaps are statistical; outer error bars without endcaps include 
systematic uncertainties. 



estimate this quantity using 

Re*- „ ^1 „ ^ (40) 

based upon Mi+ dominance and sp truncation, Figure ^2 shows that Aq T has a node near Ma and is very small for 
larger W. By contrast, A\ T peaks near Ma- Thus, it is likely that truncation errors in the multipole expansion of 
Legendre coefficients will be more serious for So+ than for Si+. This problem is illustrated in Fig. I53l which compares 
fitted values for Re(So+ /Mi+) from the multipole analysis with those based upon Eq. (1 H)[l where, because Rosenbluth 
separation is unavailable, we assume that Aq w Aq +t because Mi + dominance predicts Aq = 0. We also show 
MAID2003 calculations for both quantities. We find that MAID2003 describes both the steep slope in Re(S'o+/Mi + ) 
and the more complicated shape of 2Aq T /Aq +t fairly well, but that these quantities are rather different even in the 
immediate vicinity of W = Ma — the Legendre analysis does not even give the correct sign for this multipole ratio 
at W = 1.232 GeV. The sign difference between these quantities using MAID2003 calculations for Q 2 = 1 (GeV/c) 2 
was previously noted in Ref. |72^ and here the same problem is observed in data. The analysis of a recent experiment 
for Q 2 = 0.2 (GeV/c) 2 that measured left-right cross section asymmetries for 8 n = 20° and 160° also observed a large 
difference between ratios based upon Mi + dominance and sp truncation and those obtained by scaling MAID2003 
So+ an d Si+ multipoles to fit the data |77|. With a much more complete data set, our multipole analysis does not 
rely upon models like MAID and gives Re(S 0+ /Mi + ) = (6.4 ± 0.7)% at (W, Q 2 ) = (1.23, 1.0) directly. Assuming that 
ReMi+ w for W w Ma, the ratio Re(So+/Mi+) w ImS'o+/ImMi+ shows that ImS'o+ is positive and somewhat 
larger than the MAID2003 prediction for W = 1.23 GeV, as shown in Fig. 031 Although there is nothing special 
about M A for S 0+ , we can use the observed slope to estimate Re(S 0+ /Mi + ) = (7.1 ± 0.8)% at (W, Q 2 ) = (1.232, 1.0) 
for comparison with similar analyses purportedly at W — Ma; however, the energy dependence is steep enough that 
kinematical uncertainties could become important. Recognizing that the Q 2 dependence is mild in most models, 



imp< 

neither the present result nor that of Ref. [77| supports Schmieden's hypothesis of a large negative value for this ratio. 
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FIG. 52: (Color online) Fitted Legendre coefficients for _R£ T are compared with MAID2003 (red solid), DMT (green dashed), 
SAID (blue dash-dotted), and SL (cyan dotted). In the simplest approximation, the top panel is ReSi_ImMi+ and the second 
is Im5o + Mi + . Inner error bars with endcaps are statistical; outer error bars without endcaps include systematic uncertainties. 
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FIG. 53: Comparison between multipole and Legendre analyses of Re(5o+/Mi+). Solid and dashed curves show MAID2003 
calculations using multipole amplitudes or Legendre coefficients, respectively. Solid (open) circles show experimental ratios 
based upon multipole (Legendre) fits. The vertical dashed line indicates W — 1.232 GeV. 
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TABLE X: xt f° r Legendre and multipole analyses. 



W (GeV) xt> Legendre xt, multipole 

L17 132 L24 

1.19 1.69 1.67 

1.21 1.32 1.39 

1.23 1.50 1.65 

1.25 1.87 1.94 

1.27 1.59 1.58 

1.29 1.53 1.52 

1.31 1.61 1.42 

1.33 1.53 1.32 

1.35 1.41 1.30 



F. Interpretation of xt 

The Legendre and multipole analyses employ data for differential cross section, beam analyzing power, and recoil- 
polarization response functions with uncertainties that are primarily statistical. The cross section data include un- 
certainties in acceptance that can also be considered statistical because they are estimated from the flatness of a 
yield/simulation plateau. The uncertainties for recoil-polarization response functions are based upon diagonal el- 
ements of the covariance matrix for the maximum likelihood method. However, the fact that reduced chi-square 
values, x\,i are consistently larger than unity for both Legendre and multipole analyses suggests that uncertainties in 
extracted quantities may be underestimated. These statistics are listed in Table [X]for each W. There are several pos- 
sible explanations for this observation. First, the various recoil-polarization response functions in a given (x, W, Q 2 ) 
bin are correlated with each other, but those correlations are not considered by the Legendre or multipole analyses 
because we have no efficient means to account for them. Thus, the same fluctuation can affect several data points 
and artificially increase its contribution to \ 2 without necessarily affecting the quality of the fit. Second, systematic 
uncertainties that vary between kinematical bins were not included in the uncertainties that were used in the Leg- 
endre and multipole analyses because their effects upon various response functions are also highly correlated. Third, 
inaccuracies in baseline calculations of fixed amplitudes would impose a lower limit on x 2 even if all experimental 
correlations could be handled properly. Finally, no corrections have been made for polarized radiative corrections. 

Radiative corrections for the beam asymmetry in the p(e, e'p)w° reaction have been evaluated for Q 2 — 0.4 (GeV/c) 2 
by Afanasev et al. [f^ and found to be quite small across the A(1232) resonance. Radiative corrections for polarized 
target asymmetries are presently under investigation and generally appear to be small also |79j . but procedures for 
recoil polarization are not yet available in a form suitable for the present analysis. In principle, external radiation 
permits additional kinematical dependencies that cannot be accommodated by the response function expansions 
given in Eq. @ . Analysis of such effects probably requires an iterative procedure that begins with the current results 
to obtain model response functions, then calculates radiatively-corrected polarizations for each experimental event 
as input to an extended version of the likelihood analysis that would use a more general representation of the <f> 
dependence. In the future it may be possible to improve upon the current multipole results by iteration within a 
model of radiative corrections, and hopefully reduce xt? but that is obviously a very ambitious project. 

The simplest method of correcting for underestimates of experimental uncertainties is to multiply the uncertainties 
in extracted quantities by \[x 2 ,- We have not performed that operation here because it is somewhat arbitrary, 
assuming that neglected errors are random and uniform, but we provide Table IXl for the user's convenience. However, 
if that procedure is applied, the systematic uncertainties should probably be reduced to avoid double-counting of 
random errors presently labeled systematic. 

VIII. UNITARY ISOBAR MODEL 

Further insight can be obtained by comparing the fitted multipole amplitudes with calculations based upon the 
unitary isobar model (UIM) . We use a unitarization prescription suggested by Olssen in which the Born amplitudes 
are interpreted as contributions to the if-matrix while isobar contributions are unitarized separately using an empirical 
phase. This procedure is applied to the multipole amplitudes for isospin states. Thus, each multipole amplitude is 
expressed in the form 



A= (1 + it)B + 



(41) 
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W [GeV] W [GeV] W [GeV] W [GeV] 

FIG. 54: Comparison of fitted multipoles with UIM calculations. Inner error bars with endcaps are statistical; outer error bars 
without endcaps include systematic uncertainties. Solid curves show our UIM results while dashed curves show MAID2003. 
Dash-dotted versus dotted curves compare backgrounds for pseudovector versus pseudoscalar nNN coupling, but dash-dotted 
and solid curves for 0+ and ReMi_ amplitudes coincide. 

where B is the Born contribution (including to and p exchange), t is the ttN partial- wave amplitude, R is a Breit- 
Wigner resonance, and ip is an energy-dependent phase. Angular momentum and isospin labels have been suppressed 
for brevity. The unitarity phase is adjusted according to the Fermi- Watson theorem |65|. which requires 

Be** =±\R\e u (42) 

where S is the ttN elastic phase shift and is assumed to be real throughout our energy region. These phases differ 
for each multipole and depend upon the background parametrization. The sign ambiguity is resolved locally and is 
significant only for the P\\ partial wave at low energies. 

We used the SAID FA02 solutions for ttN partial- wave amplitudes Resonances were parametrized in Breit- 
Wigner form, Eq. (|3 1 11 - with energy-dependent widths based upon the MAID model Only the A and Roper 

resonances are significant in our W range. Widths and branching ratios were taken from the most recent compilation 
from the Particle Data Group from which widths at resonance were taken as 120 MeV for the A and 350 MeV, 
with a ttN branch of 65%, for the Roper. 

UIM calculations are compared with fitted multipole amplitudes in Fig. [SU The dash-dotted versus dotted curves 
compare background amplitudes for pseudovector (PV) versus pseudoscalar (PS) nNN coupling and differ only for 
0+ and 1— multipoles. The solid curves show the total UIM multipoles for pseudovector coupling, which coincide with 
dash-dotted curves for 0+ and ReMi_ amplitudes because no resonances were included for those amplitudes. The 
dashed curves show MAID2003, which uses a similar unitarization prescription but interpolates between pseudovector 
coupling for small W and pseudoscalar coupling for large W |67| . The background amplitudes acquire their imaginary 
parts from the irN phase shift. This effect is especially important for the 1+ and 0+ multipoles and is responsible for 
the differences between the real parts of background amplitudes shown in Fig. 1 5 41 and the baseline amplitudes shown in 
Figs. 1421451 The effect of the -P33 resonance upon M\ + Born amplitudes is particularly strong. The best sensitivity to 
the wNN coupling is in real parts of 0+ multipoles, especially ReSo+. The data for ReS'o+ clearly favor pseudovector 
over pseudoscalar coupling; the admixture used by MAID2003 flattens the W dependence but is not supported by 
these data. The data for ReEo + also prefer pseudovector coupling but could accommodate the MAID2003 admixture 
if a sufficiently strong 511(1535) contribution were included; we did not investigate that possibility because our data 
are limited to relatively small W. 

Neither our version of the UIM nor that of MAID2003 reproduces the imaginary parts of 0+ amplitudes well. Both 
approximate the average value for Im-Eo+ but neither reproduces its W dependence. Note that Imi?o+ is especially 
sensitive to the axial form factor through the unitarization factor. We used a dipole form factor with Ma = 1.0 
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(GeV/c) 2 . The discrepancy is more severe for ImS'o+ which grows much more rapidly than either UIM, at least for 
PV coupling. One might be tempted to include a strong Su contribution and to adjust the PS/PV admixture to 
maintain the fit to Re5o+, but we are loath to attempt such a fit without data across, or at least closer to, the 
5ii(1535) resonance. Therefore, our full calculations, shown as solid lines, assume pure pseudovector ttNN coupling 
and do not include an Su resonance. The discrepancy in ImS'o+ bears further investigation. 

The present UIM reproduces ImMi+ very well, but appears to be missing a small positive contribution to ReMi+. 
Nevertheless, the accuracy is sufficient to quote a resonance contribution of M^ 2) = (2.91 ± 0.15) p.b 1 ' 2 where 
the estimated uncertainty is qualitative. The E\+ and Si+ calculations use the quadrupole ratios from Sec. IVI Dl 
without further adjustment. These calculations are qualitatively consistent with the data but do not reproduce 
the W dependence accurately. Interestingly, the present UIM is more successful for lmEi + while MAID2003 is more 
successful for ImS'i + . It might be possible to improve these fits by adding phenomenological background contributions 
to the B terms of Eq. lj¥T|) . 

Although our W acceptance does not span the Roper resonance, its width is broad enough to permit estimation of 
the electromagnetic helicity amplitudes based upon data for its low-W side. Having selected pure PV coupling, we 
find little room in ReAfi_ for transverse excitation of the Roper resonance, such that p A 1 / 2 is consistent with zero at 
Q 2 = I (GeV/c) 2 . That conclusion obviously depends upon the PS/PV admixture, with more PS coupling permitting 
larger p A 1 / 2 . Conversely, reproduction of Re<Si_ requires appreciable longitudinal Roper excitation, larger for PV 

than for PS coupling. With PV coupling, we estimate p Si/ 2 = (0.05 ± 0.01) GeV -1 / 2 where the uncertainty is again 
qualitative, based upon simultaneous fits to both real and imaginary parts. Similarly, Laveissiere et al. |2Cj | fit cross 
section data for Q 2 — 1.0 (GeV/c) 2 by reducing the MAID2003 estimate of p A 1 / 2 to a value consistent with zero but 
their estimate of p 5i/ 2 = 0.019±0.010 GeV~ 1/2 is smaller than ours; however, it should be noted that their angular 
range was quite limited. These findings are also qualitatively consistent with those of Aznauryan et al. 82[ based 
upon cross section and R' LT data for the pn and nir + channels at Q 2 = 0.4 and 0.65 and pr\ cross section data for 
Q 2 = 0.375 and 0.75 (GeV/c) 2 with W spanning the second resonance region. They found that p A 1 / 2 is small and 
appears to change sign near 0.5 (GeV/c) 2 . Our value for p S\/ 2 is somewhat larger than theirs but some models, such 
as those by Capstick and Keister [8j|, feature a peak in p S\/ 2 for Q 2 ~ 0.7. Furthermore, unlike the data set employed 
by Ref . |82| , recoil polarization provides sufficient phase information to separate multipolc amplitudes explicitly. It is 
obviously desirable to acquire new data for either recoil or target polarization that reach larger W. Nevertheless, the 
current results and those of Ref. appear to exclude the hybrid baryon model of the Roper resonance for which 
electromagnetic excitation would be purely transverse 0, 0] • 

Clearly it would be of interest to perform an energy-dependent analysis of the entire data simultaneously that uses 
irN phase shifts to enforce unitarity. However, the unitarization procedure is not unique [84[ and requires background 
models for both pir° and nir + . Furthermore, Born diagrams at tree level do not necessarily represent background 
amplitudes with sufficient accuracy. An even more ambitious analysis could use dispersion relations to improve the 
background model. However, such an analysis should consider all available data and is beyond the scope of the present 
work. Therefore, we have not attempted to optimize the UIM parameters, except for A and Roper strengths. 

The present multipole analysis can be described as energy-independent because each W is fit independently. It also 
uses the minimum possible theoretical information, simply the small Born amplitudes for partial waves with > 2. 
No attempt has been made to enforce unitarity in the multipole analysis and some concern has been expressed that 
the steep slope we find for ImSo+ might be inconsistent with the unitary isobar model [slf. On the other, the good 
agreement between SAID calculations and the data for both real and imaginary parts of 5*0+ shows that it is possible 
to describe the pir° data for this multipole in a manner that is consistent with unitarity (see Fig. I43f) . An experimental 
test of unitarity that is independent of models for resonant and nonresonant amplitudes would require comparable 
polarization data for the n7r + reaction, but such data are presently unavailable. 

IX. SUMMARY AND CONCLUSIONS 

We measured angular distributions for differential cross section, beam analyzing power, and recoil polarization in 
the p(e,e'p)TT° reaction at Q 2 = 1 (GeV/c) 2 with 1.17 < W < 1.35 GeV across the A resonance and have obtained 
14 separated response functions and 2 Rosenbluth combinations, of which 12 have been measured for the first time. 

We compared the data for response functions with calculations for four recent models: MAID, DMT, SAID, and 
SL. Variations among these models are relatively small at W ss Ma for quantities that depend upon real parts of 
interference products, but increase with W. Variations among models are much larger for quantities dependent upon 
imaginary parts that are more sensitive to background amplitudes. MAID and DMT are similar and in relatively 
good agreement with data for W ~ Ma, but neither provides a uniformly good description of the data for larger W. 
The SL model, which does not include higher resonances, underpredicts the cross section for larger W while DMT 
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is too strong. The SAID model has considerable difficulty with helicity-independent LT response functions that are 
probably caused mostly by its rather strong ReSi_ amplitude. 

We performed a multipole analysis that fits both real and imaginary parts of the multipole amplitudes for low 
partial waves while those for higher partial waves are constrained by either Born terms or by the best available model 
calculations. Fitted multipole amplitudes based upon Born, MAID, DMT, or SL models arc practically indistinguish- 
able, but the available version of SAID does not provide a suitable baseline because some of its > 2 amplitudes 
are too strong. The final analysis is based upon the Born model to minimize bias. We chose not to vary Im_Mi_ 
in the final analysis because it is predicted to be negligible in our energy range but its fitted values are strongly 
correlated with those of Im5i_ for the present data set. We were able to extract consistent results for all 1^ < 1 
amplitudes, except ImMi_, plus the real parts of 2— multipoles. The most significant differences between fitted and 
model amplitudes arc found in 0+ and 1— multipoles. The data also show that ImS'o+ grows faster than predicted by 
MAID, DMT, or SL, but is described reasonably well by SAID. Good sensitivity to ReMi_ is provided by the Rfp T 
response function; there is a wide spread among models and MAID2003 fits the ReMi_ data best and is close to the 
Born baseline. Similarly, the best sensitivity to ReSi_ is provided by i?£ T but none of the models is accurate — SAID 
is much too strong while MAID, DMT, and SL are too weak for that amplitude. The data are substantially stronger 
than the Born amplitude, suggesting significant longitudinal Roper contributions arising from a radial excitation. 

We find that truncation errors in the traditional Legendre analysis of N — > A quadrupole ratios can be significantly 

larger than statistical errors. Using parabolic fits to the energy dependence, we obtain RgZr ' = (— 5.87± 0.20)% and 

Rem ] = (-1-76±0.19)% from the traditional analysis or i?^ 0) = (-6.61±0.18)% and R^ ] = (-2.87± 0.19)% from 
the multipole analysis for W = 1.232 GeV and Q 2 = 1.0 (GeV/c) 2 . These results are consistent with the single-energy 
analysis published previously [33|. The model dependence of the multipole analysis is small and the Legendre fits 
are stable with respect to the number of fitted terms, yet the differences between these analyses are several standard 
deviations. We have demonstrated that the multipole analysis is more reliable because it does not depend upon M\ + 
dominance or sp truncation. Both model calculations and the multipole analysis of data demonstrate that neither 
assumption is reliable and that multipole products omitted by that truncation scheme make important contributions 
to the Legendre coefficients that spoil the accuracy of the simple estimators of quadrupole ratios employed by the 
traditional Legendre analysis. Truncation errors are especially severe for Rem- 

We also find that Re(5o+/Mi-|_) = (7.1 ± 0.8)% at W = 1.232 GeV is qualitatively consistent with most recent 
models and with a recent measurement j?3 at Q 2 — 0.2 (GeV/c) 2 of left-right cross section asymmetries at a pair of 
supplementary proton angles, but is inconsistent with a recent hypothesis |22| that a large negative value is needed to 
explain inconsistencies between SMR analyses at Q 2 = 0.13 (GeV/c) 2 using earlier data for forward versus backward 
On- Truncation errors in the Legendre estimator for Re(So+/Mi + ) are quite severe [72, even resulting in an incorrect 
sign at Q 2 ~ 1 (GeV/c) 2 . The analysis in Ref. jf?} relied on the MAID model instead of the Legendre estimator but 
accurate, model-independent results require a phase-sensitive multipole analysis as performed here. 

Finally, we compared the fitted multipole amplitudes with calculations based on a unitary isobar model (UIM). 
The ReSo+ and Rei?o+ multipoles strongly prefer pseudovector over pseudoscalar irNN coupling and do not support 
the proposed mixing between PV and PS coupling employed by the MAID model. However, the UIM does not 
reproduce imaginary parts of 0+ multipoles well, with Im5o+ increasing more rapidly with W than expected from 
the pseudovector Born contribution. Note that the ability of SAID to reproduce both real and imaginary 5*0+ data 
for pir° shows that failure of UIM to reproduce these data does not necessarily mean inconsistency with unitarity; a 
rigorous test of unitarity would require similar nir + data. The UIM reproduces M\ + well, but the W dependencies 
of Imi?i + and ImSi+ are only qualitatively consistent with the data. If we select pure PV coupling, the data for 
ReMi_ are described well by Born terms, suggesting that p Ai/ 2 for the Roper is consistent with zero. Therefore, 
Roper excitation is dominantly longitudinal for Q 2 = 1.0 (GeV/c) 2 , where we find p S\/2 = (0.05 ± 0.01) GeV~ 1//2 . 
Although a larger PS admixture would permit appreciable transverse Roper excitation, these findings tend to exclude 
the hybrid baryon model of the Roper resonance. Clearly it would be of interest to extend these measurements to 
larger W. 

In conclusion, recoil and/or target polarization data are essential to multipole analyses of meson electroproduction 
reactions, providing access to the relative phase between resonant and nonresonant contributions. Although neutral 
pion electroproduction in the A region is the easiest example, this experiment demonstrates the feasibility of the 
method and we hope that it will be applied over wider kinematic ranges and to related reactions. An advantage of 
this type of analysis is that it minimizes the dependence upon models; however, it does not guarantee that the fitted 
multipole amplitudes will depend smoothly on both W and Q 2 . Model-dependent analyses which adjust parameters 
of an effective Lagrangian or unitary isobar model should produce kinematically smooth multipole amplitudes at the 
expense of possible bias. Presumably, analyses of these types would also be less sensitive to variations of acceptance- 
averaged W and Q 2 between bins of the angular variables, (x,(j>). Both types of analyses would benefit from more 
extensive coverage in W. With sufficient kinematic coverage one hopes to obtain reliable transition form factors for 
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overlapping resonances. 
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APPENDIX A: SPIN TRANSPORT 

The spin transport matrix consists of the sequence of transformations 

F = Rfpp C -Rspectrometer -Shall R\V T = ST (Al) 

where 

• Rw performs an active Wigner rotation of the polarization vector from the center of mass to the laboratory 
frame, 

• -Rhaii performs a passive transformation from the reaction basis to the hall basis, 

• ^spectrometer transforms from the hall basis to the spectrometer basis, 

• C transports the spin through the magnetic spectrometer, and 

• i?fpp transforms to the local FPP coordinate system. 
The individual transformations are detailed below. 



1. Wigner rotation 



The polarization vector is transformed from the cm reaction basis with £ along p cm to the laboratory frame with L 
along plab using a rotation of the form [8?" 

cos 9w sin 9w 

1 I I P n I (A2) 




— sin 6w cos 9 




where 

taxi0 w = — — - (A3) 

7cm (P COSy 

cm T" Pcm ) 

Here (3 cm = p C m/E cm is the nucleon velocity in the cm and (3 = q/ (m p + u) is the velocity of the cm relative to the 
lab. This matrix is identified as Rw in Eq. IjAlfl . 



2. Transformation to spectrometer frame 



The hall basis is defined with z along the beam line and y vertically upward in the lab. It is useful here to define 
the horizontal angle a in the xz plane measured counterclockwise from z and the vertical angle (3 to be positive above 
the horizontal plane. Unit vectors q along the momentum transfer and p along the nucleon momentum plane then 
take the form 

q = (sin a q cos (3 q , sin ft q , cos a q cos (3 q ) ( A4a) 

p = (sina p cos/3 p , sin/3 p , cosap cos/3 p ) (A4b) 
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where both a p and a q are negative for our configuration with the proton spectrometer on the right side of the beam. 
The laboratory reaction basis (S, N, L) is now defined by 

L = p (A5a) 

N oc qAp (A5b) 

S = N A L (A5c) 

where L is along the nucleon momentum, N is normal to the reaction plane, and S is sideways within that plane. 
Therefore, the transformation from the laboratory reaction to the hall frame employs 



flhaii = [S, N, Lj (A6) 

obtained by usin g th ese unit vectors as the columns of a square matrix. 

The COSY [871 l88j calculation of spin precession employs the TRANSPORT coordinate system with z s along 
the central axis of the spectrometer, x s downward, and y s toward the left of the spectrometer midplane. Thus, 
transformation from the hall to the spectrometer frames is accomplished using 

-1 

-^spectrometer = I COS « - sin O> I (A7) 




sin O.Q 

where the central angle of the spectrometer, ao, is also negative for our configuration. 

3. Spin transport 

The spin transport matrix dj relates components of spin in the final spectrometer basis to those in the initial 
spectrometer basis, where these bases differ by a rotation about the y s axis through the bend angle Hq for the central 
ray. Thus, the spin transport matrix C^J for a pure dipole would take the form 

/ cosxo -sinxo \ 
C (0) =01 (A8) 
y sin xo cos xo J 

where xo — JK-p^o, k p is the anomalous magnetic moment, and 7 = E/m. The central bend angle is tto = 45° for 
HRS. 

Elements of the spin-transport matrix were computed by the differential-algebra transport code COSY using a 
magnetic model of HRS and were expanded in the form 

C ij= C^ np x k 9 e y m 4> n 6 p K (A9) 

k,l,m,n,p 

where (x,0,y,<fi) are the reconstructed track variables expressed in TRANSPORT form and 8k = (K — K )/K is the 
kinetic energy displacement relative to the central value. The expansions are carried to 5 th order. Further details can 
be found in Ref. 

4. Transformation to FPP basis 

Tracks before and after scattering by the analyzer are reconstructed in the form 

ki — (sina^ cos/3j, sin/^i, cos cos/?i) (AlOa) 
kj = (sin a/ cos /?/, sin/3/, cosa/ cos/?/) (AlOb) 

where a and (3 are cartesian angles. The incident polarization is transformed to a basis aligned with ki using 

cos on — sin <Xi 

Rfpp = \ — sin a>i sin fa cos fa — cos ctj sin fa | (All) 
sin at cos fa sin fa cos a« cos fa 
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Similarly, the polar and azimuthal scattering angles in the FPP are determined using 

sin fpp cos </>fpp \ 

sin 6>f pp sin 0f pp \ = R ipp k f (A12) 
cos 0{ pp J 

such that 

F ■ n fpp = F 2 cos 0f pp - Fi sin fpp (A13) 
where nf pp = h A kf /\ki f\ kf \ is a, unit vector normal to the FPP scattering plane. 

APPENDIX B: MAXIMIZATION OF LIKELIHOOD 

The likelihood function for polarization measurements takes the generic form 

L = H^-(a i + X i -R) (Bl) 

i 

where the Latin index i enumerates events satisfying the selection criteria for a particular kinematical bin, R is an 
n-dimensional result vector where R a represents an acceptance-averaged quantity, identified by the Greek index a, 
that is assumed to be constant for all events in the bin, X ia is an event-dependent coefficient that weights the effect of 
parameter R a , and m includes the effects of instrumental asymmetry and is assumed to be of order unity. The scalar 
product denotes contraction with respect to a. The logarithm of the likelihood function is maximized with respect 
to R a when 

y = 

i 

Although this set of n equations cannot be solved in closed form, an iterative solution can be obtained by combining 
successive substitution with linearization of the summand. 

If the asymmetry is small, linearization provides an approximate solution in the form 

R w A -1 • V (B3) 

where 

K. = E7 ( B4 ) 

u>i 



is a measurement vector and 



a; a. 



is the design matrix for the experiment. More generally, let 

R (k) = R (k-i) + AR (k) ^ 

represent a sequence of improved approximations to the result vector, and let 

represent the measurement vector and design matrix for iteration n. Thus, we obtain 

Ai?W = (A**)) -1 -V^ (B8) 
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When the asymmetry is small, one normally begins with the unbiased initial estimates Ra^ = and expects rapid 
convergence, but if the asymmetry is large it may be advantageous to begin with model estimates for the 
parameters because each iteration for a large sample may require considerable computation time. 

Let R represent the maximum-likelihood estimates of the model parameters given by the convergence of Eqs. (|B6I - 
IB8|) and let AR = R— R represent the deviation vector. We assume that the likelihood function is described well by 
the Gaussian approximation 

lnLwlogio-iAi?- A- AR (B9) 



near its maximum. Therefore, the covariance matrix is given by a — A 1 where 

(a, + Xt ■ R) 2 



= £ . (BIO) 



and we estimate the parameter uncertainties 



as standard deviations. 



on a = o)iZ (Bll) 



APPENDIX C: AMPLITUDES 

The reaction amplitudes for any A(e,e'N)B process where A has spin-i and B spin-0 that is governed by the 
one-photon exchange mechanism can be expressed in terms of helicity amplitudes of the form 

nx f x^(W,Q 2 ,9^) = (A/I^IA*,^) (CI) 

where Aj and A/ are the initial and final helicities of the nucleoli, A 7 is the helicity of the virtual photon, is an 
appropriately normalized transition current operator, and e M is the virtual-photon polarization vector. Since parity 
conservation [gj| requires \H.—\,—\ i —\A = |Wa/A 4 A 7 |) it is sufficient to consider six independent amplitudes Hi for 
(Ajj , Aj, A 7 ) chosen as (— — \, 1), (— \, \ , 1), — , 1), (|, \ , 1), {\, \ , 0), and (|, — |, 0) and numbered sequentially 
90, 9]J. The most general current operator for pseudoscalar meson production consistent with parity conservation 
and other symmetries can be represented in term of CGLN amplitudes 92, 93] as 



iT a = -{F' b (T ■ q + F' & (j -p) (C2a) 

UJ 

IT = T\S — iF<xa ■ pa x q + T 3 pa ■ q + T^pS ■ p + T^qa ■ q + T^qa ■ p (C2b) 

where 

T' b = H+T z p-q + T x (C3a) 
T' 6 = Ts + T A p-q (C3b) 

are combinations which simplify the multipole analysis. Using phases for helicity states following the conventions of 
Jacob and Wick [g^, the helicity amplitudes are related to CGLN coefficients by 

Q 

Hi = -e l0 v / 2sin6»cos(-)(.7 : 3 + T A ) (C4a) 
H 2 = V2 ^cos(^)(^ 2 -^ 1 )+sin0sin(^)(^ 3 -^4)^) (C4b) 

H 3 = e 2 ^-±= sin 9 sin (t) [T z -Fi) (C4c) 

/ 9 9 \ 

Hi = x/2(sm(-)(F 1 +F 2 )+sm9cos(-)(F 3 +F i )j (C4d) 

H 5 = ^cos(^)(^+^) (C4e) 

UJ z 

n 6 = e ^ sin A ( ^_^ 6 ) (C 4f) 

UJ L 

(C4g) 
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where 9 and <j) refer to the meson. 

CGLN amplitudes can be expanded in terms of multipole amplitudes as follows [9^ |. 

T x = ^{tM t+ + Ei + )Pl, +1 {x) + {{i + l)M t -+E t -)P' l _ 1 {x) 
i 

i 

i 

T 4 = J2( M *+ - E *+ - M t- - E e -) P' t '{x) 
i 

J2(tSe--(l + l)St + )Pl(x) 



T' 

•J R 



U) 

1 



The multipole amplitudes can be projected from CGLN amplitudes using [9C 



M l+ 
M t - 
Ee+ 
Eg- 
Si+ 
Si- 



1 



dx 



1 



2£ + 1 



{P&+X — Pe-i)J~3 



= h 1 ,Lr 



—P(Ti + Pg-\T<2, — 



1 



2£ + 1 



2(£+l) J-x 



dx 



PiTx - Pi+X^2 



(Pe+x — P(.-xh3 



(P e+1 -Pe-x)F 3 



2£+l 



1+1 
2^ + 3 



{P t - P e+2 )F 4 



hL dx 

q 1 

w2(i+i) y_ 



1±1 

21+1 



21- 1 



2f 



y [P^ + Pi+xFq] 

J dx [P t F' 5 + Pt-iFt] 



(C5a) 
(C5b) 
(C5c) 
(C5d) 
(C5e) 
(C5f) 

(C6a) 
(C6b) 
(C6c) 
(C6d) 
(C6e) 
(C6f) 



where x — cos 0^ • 

Often it is useful to express multipole amplitudes in terms of partial-wave helicity amplitudes |90l | according to 



Mm 
£m 
£m 
M e . 
Et- 



where 



[2A t+ -(£+ 2)B l+ ] I [2(1 + 1)] 
[2A e+ +£B I+ ]/[2(£ + l)] 
Ct+/{1+1) 

[2A e - + (l-l)B e -}/(2£) 
[-2At- + [£ + l)B l J\/(2£) 
-Ct-ll 



~ [*Af /+ + (i + 2)P,+] 



P^+ 



£!^_i_ — Mi 
(i + l)S e+ 



1 



[(*+l)Af<_ -(£-l)E e _] 



E t - + M t - 
-IS.- 



(C7a) 
(C7b) 
(C7c) 
(C7d) 
(C7e) 
(C7f) 



(C8a) 

(C8b) 
(C8c) 

(C8d) 

(C8e) 
(C8f) 



differ from the usual Hebb- Walker convention |9l| by using Si± in place of L(±. 
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APPENDIX D: RESPONSE FUNCTIONS 

The response functions can be expressed in terms of helicity amplitudes Hi(W, Q 2 ,0) = Hi(W, Q 2 , 8, 0) for <j> = 
as follows. 

R l = f!> (W + lffel 2 ) (Dla) 
Q 2 

sm0 Rl = -2j-lmH 5 H£ (Dlb) 

Rt = l -{\H l \ 2 + \H 2 \ 2 + \H 3 \ 2 + \H 4 \ 2 ) (Die) 

sm6R$ = 1111(^^+^4*) (Did) 

sin0i? LT = -j^Rc^ -H 4 )H*-(H 2 +H 3 )H*) (Die) 



i?2 T = -^ Im «^ +(H 1 -H 4 )H*) (DID 



sin i?i T = -^Im ((ifx + tf 4 )ff 5 * - (H 2 - H 3 )H* 6 ) (Dig) 

R\ T = -^Im ((tf 2 - ff 3 )^5 + (#i + ^4)^ 6 *) (Dlh) 

sinfli?^ = --^Im((H 1 -H 4 )^*-( J ff 2 + J ff 3 )i/*) (Dli) 

R'E T = -^^ c ((^2 + H 3 )H* + (H, - H 4 )H* 6 ) (Dlj) 

sm0R'l T = --X-R e ((ff 1+ ff 4 )ff* -(H 2 -H 3 )H£) (Dlk) 

i?2 T = --^Re((ff 2 -H3)ff5* + (^i+^4)H 6 *) (Dll) 

sin 2 6 R TT = Re (H 2 H 3 — HiH 4 ) (Dim) 

sin(9i?£ T = Im(i?iF 2 * + H 3 H* 4 ) (Din) 

sin 2 6» i?^T = Im (#2#3 - H^l) (Dlo) 

sin<9i^ T = lm(H 3 Hl- H^H*) (Dip) 

R'St = ^ (l#3| 2 + |^4| 2 - l^ll 2 - |^ 2 | 2 ) (Dlq) 

sin 6>i?£ T = Re(i2i#3 + i? 2 tf 4 ) (Dlr) 

APPENDIX E: BORN BASELINE MODEL 

In this section we summarize the Born baseline model used for the multipole analysis, including only the terms that 
contribute to the pn° channel. The electromagnetic vertices are represented by effective lagrangians of the form 



-7AW 



= -eN 



N (Ela) 



C^v = e— e^dM^^fe^ 3 + /?f)^v(Q 2 ) (Elb) 

(Elc) 

where N represents a nucleon field operator, is the electromagnetic vector potential, 7r is the pion field as an 
isospin vector, V G W,p} denotes a vector meson. We used conventional dipole and Galster form factors for the 
nucleon and monopole form factors 

F,MQ 2 )= ( 1 + ^Pj (E2) 
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TABLE XI: Parameters for vector-meson vertices. 



V 


m v [MeV] 


gv 1 


gv 2 


Avnn [MeV] 


Xv 




782.6 


21 


-12 


1200 


0.314 


P 


769.0 


2 


13 


1500 


0.103 



for jirV vertices. The jnV and VNN parameters are listed in Table IXTI and were taken from Ref. |67j . 
We used pure pseudovector nNN coupling 

C*NN = -™^757mT • (^tt)N (E3) 
zmjv 

with ff^jvjv = 13.4. Drechsel et al. [(57j proposed a more flexible ttNN model that interpolates between pseudovector 
coupling for small p v and pseudoscalar coupling for large p v , but this variation only affects real parts of 0+ and 
1— baseline multipoles and the fitted parameters simply compensate for variations of the pseudoscalar/pseudovector 
mixture anyway. 

Finally, the VNN coupling is described by 



£-VNN — —N 



NF VNN (t) (E4) 



where u> v and p v represent u) and p fields. A strong form factor, 

A 2 _™2 



F VNN (t) = m v (E5 ) 

IV VNN 1 



is applied to the VNN vertex according to the prescription of Brown et al. j^H . 
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